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NUMERICAL  MODELING  OF  SILICON  PHOTODIODES 
FOR  HIGH-ACCURACY  APPLICATIONS 

Introduction  to  NISTIR  4592 


The  purpose  of  this  National  Institute  of  Standards  and  Technology  Interagency/Internal 
Report  (NISTIR)  is  to  provide  the  source  code  for  two  Turbo  Pascal*  5.5  programs  and  an 
MSDOS  batch  program,  along  with  a paper  [l]  that  describes  the  programs  and  provides 
examples  of  their  use.  These  three  programs  serve  as  a batch  mode  interface  to  support 
high-accuracy  photodiode  modeling  with  Version  2 of  the  semiconductor  device  modeling 
program  PC-lD  [2].  These  programs  are  useful  because  the  interactive  user  interface  of 
PC-lD  is  optimized  for  solar  cell  modeling,  and  it  is  somewhat  difficult  to  access  the  highest 
accuracy  available  from  PC-lD  through  this  interface. 

The  paper  included  in  this  NISTIR  is  divided  into  Parts  I,  II,  and  III.  Part  I describes  PC- 
lD  from  the  point  of  view  of  high-accuracy  photodiode  modeling,  and  it  also  describes  the 
programs  that  support  its  use  in  this  application.  The  source  code  listed  in  the  appendices 
to  Part  I could  be  easily  modified  to  support  other  specialized  applications  of  PC-lD.  Parts 
II  and  III  present  examples  of  the  use  of  the  programs  described  in  Part  I to  model  two 
different  types  of  experiments  performed  on  silicon  photodiodes  in  various  high-accuracy 
applications. 
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the  particular  purposes. 
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Part  I:  SIMULATION  PROGRAMS 

Jon  Geist  and  Deane  Chandler-Horowitz 
NIST,  Gaithersburg,  MD  20899 
A.  M.  Robinson  and  C.  R.  James 
University  of  Alberta,  Alberta,  Canada  T6G  2G7 

Abstract 

The  suitability  of  the  semiconductor-device  modeling  program  PC-lD  for  high-accuracy 
simulation  of  silicon  photodiodes  is  discussed.  A set  of  user  interface  programs  optimized  to 
support  high-accuracy  batch-mode  operation  of  PC-lD  for  modeling  the  internal  quantum 
efficiency  of  photodiodes  is  also  described.  The  optimization  includes  correction  for  the 
dark  current  under  reverse-  and  forward-bias  conditions  before  calculating  the  quantum 
efficiency,  and  easy  access  to  the  highest  numerical  accuracy  available  from  PC-lD,  neither 
of  which  are  conveniently  available  with  PC-lD’s  standard  user  interface. 

Key  words:  high  accuracy;  internal  quantum  efficiency;  PC-lD;  photodiode  modeling; 
silicon  photodiodes. 


1.  Introduction 

PC-lD  [1-2]  is  a computer  program  for  numerical  modeling  of  the  electrical  performance 
of  one-dimensional  semiconductor  devices  including  photogeneration.  It  is  optimized  for 
solar  cell  modeling,  and  runs  on  IBM  PC-compatible  computers  equipped  with  a numerical 
coprocessor.  Version  2 of  PC-lD  includes  realistic  models  of  most  of  the  semiconductor 
material  and  device  properties  that  are  important  for  high-accuracy  applications  of  pho- 
todiodes. Consequently,  Version  2 of  PC-lD  has  the  potential  to  be  a tool  for  supporting 
these  applications. 

Unfortunately,  for  reasons  that  have  to  do  with  the  difference  between  the  applications  for 
high-accuracy  photodiodes  and  those  for  solar  cells,  it  is  not  easy  to  achieve  the  highest 
accuracy  of  PC-lD  through  its  standard  user  interface.  Luckily,  the  designer  of  Version  2 
of  PC-lD  anticipated  this  sort  of  problem  and  included  the  option  to  run  it  in  batch  mode 
in  a way  that  does  allow  access  to  the  highest  levels  of  numerical  precision  and  accuracy 
that  are  available  from  PC-lD. 

This  paper  reports  the  development  of  a program  shell  for  PC-lD  that  provides  a batch- 
mode user  interface  optimized  for  high-accuracy  modeling  of  photodiodes.  The  shell  con- 
sists of  three  programs.  The  first  is  a program  that  prepares  the  input  for  PC-lD,  the 
second  is  a program  that  reads  the  output  from  PC-lD,  and  the  third  is  a MS-DOS  batch 
file  program  that  supervises  the  sequential  execution  of  PC-lD  and  the  other  two  programs. 

The  remainder  of  Part  I of  this  paper  reviews  PC-lD  and  describes  the  new  program  shell 
for  high-accuracy  modeling  of  photodiodes.  Typical  applications  are  described  in  Parts  II 
and  III. 
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2.  Description  of  PC-lD 


PC-lD  is  an  interactive,  graphical,  semiconductor-device  simulation  program  that  solves 
the  fully-coupled,  drift-diffusion  (electron  and  hole  transport)  equations  in  one  dimension. 
Its  principal  limitations  are  that  it  supports  no  more  than  three  regions  of  possibly  different 
materials,  and  that  it  is  a one-dimensional  model  with  a maximum  of  150  finite  elements. 
Within  these  constraints,  the  program  is  very  versatile.  It  is  well  beyond  the  scope  of  this 
paper  to  describe  its  various  capabilities. 

Both  the  IBM  Pascal  source  code  and  executable  code  are  available.  The  latter  requires 
an  IBM  PC-compatible  computer  equipped  with  at  least  512  Kbytes  of  random  access 
memory  (RAM),  a 80X87  coprocessor  chip,  a CGA,  VGA,  or  EGA  compatible  graphics 
adapter,  and  a matching  display.  To  use  the  shell  described  in  this  paper,  batch-mode 
operation  is  necessary  which  requires  an  MS-DOS  compatible  operating  system  (Version 

3.1  or  greater). 

Only  one  of  the  interactive  modes  of  PC-lD  is  of  interest  for  the  purposes  of  this  paper, 
the  one  that  allows  the  creation  and  saving  of  parameter  files.  Parameter  files  contain 
all  of  the  information  needed  by  PC-lD  to  model  a specific,  user-defined  device.  These 
files  are  stored  in  binary  format  under  user-defined  names  with  .PRM  as  their  extension. 
Included  in  the  parameter  file  are  the  names  of  other  data  files  needed  by  PC-lD  to  model 
the  device. 

2.1  Opening  Menu 

Figure  1 shows  the  opening  menu  when  PC-lD  is  run  in  the  interactive  mode.  An  overview 
of  the  options  mentioned  in  this  menu  is  appropriate  for  what  follows.  The  first,  second,  and 
fourth  options,  “Proceed  with  Solution,”  “Preview  / Examine,”  and  “Output  Graphs:”  are 
used  in  the  interactive  modeling  mode,  and  are  of  no  concern  here.  The  “Solution  Mode:,” 
“Base  Voltage:,”  “Collector  Voltage:,”  and  “Light:”  options  are  used  to  set  the  electrical 
and  optical  conditions  under  which  the  simulation  will  be  run.  For  the  purposes  of  this 
paper,  it  does  not  matter  what  selections  are  made  for  these  options,  because  they  will  be 
set  to  the  correct  values  by  the  high-accuracy  photodiode  modeling  shell  when  it  is  run. 

The  “Simulation  Parameters:”  option  is  used  to  load  parameter  files  in  the  interactive 
mode,  but  it  is  also  used  to  save  parameter  files,  which  is  of  interest  to  the  use  of  PC- 
lD  in  batch  mode.  The  “Reinitialize:”  and  “Number  of  Finite  Elements:”  options  allow 
details  of  the  numerical  procedures  to  be  modified  to  a certain  extent.  Normally,  150 
finite  elements,  the  maximum  allowed,  will  be  needed  for  high-accuracy  modeling.  It 
will  sometimes  be  necessary  to  use  this  menu  option  in  order  to  force  PC-lD  to  use  all 
150  elements.  Reinitialization  is  probably  desirable  except,  perhaps,  when  simulating  a 
reverse-bias  experiment. 

The  remaining  options  define  the  detailed  nature  of  the  device  and  external  circuit  to  be 
modeled.  All  of  the  parameters  obtained  from  the  internal  models  built  into  PC-lD  are 
adjusted  to  the  values  appropriate  for  the  temperature  entered  following  the  “Tempera- 
ture:” heading  at  the  bottom  of  the  opening  menu.  The  user,  of  course,  must  assure  that 
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PC-ID  PROGRAM  CONTROL 


l}JW 

^view  / Imm  (None) 
hnulaiion  ParaMe ten's;  1337 JEU 
Output  Graphs;  Specified 
Reinitialize;  Ves 
NuMber  of  Finite  Elenents;  150 
Solution  Mode;  EquilibriuH 
TeMperature;  25  C 
Area;  Uni  font  1 cm2 
Thicknesses;  1 un,  299  un 
Materials;  Sl«,  SI* 

Doping;  1337,  Internal 
RecoMbination;  Internal 
Surfaces;  Charged,  Neutral 
Circuit;  Ubb 
Base  Voltage;  Zero 
Collector  Voltage;  Zero 
Light;  10  uN^mR 

Photogeneration;  MonochroHe  440  nM 
Reflectance;  0,  0.85 


Figure  1.  The  opening  menu  of  PC-lD  when  run  in  the  interactive  mode. 


any  user-defined  parameters  or  data  files  are  consistent  with  the  chosen  temperature.  The 
“Area:”  option  is  straightforward.  The  “Thickness:”  option  allows  three  different  regions 
to  be  defined,  each  having  its  own  thickness,  material,  and  doping  properties.  If  the  thick- 
ness of  a region  is  set  to  zero,  it  is  not  included  in  the  solution,  and  it  is  not  explicitly 
indicated  following  the  “Thickness:”  heading. 

2.1.1  Materials  Properties 

Selection  of  the  “Materials:”  heading  of  the  opening  menu  allows  such  properties  as  the 
carrier  mobilities,  band-gap  narrowing.  Auger  recombination  cross  sections,  and  absorption 
coefficient  to  be  adjusted.  For  each  region,  a set  of  three  files  having  the  extensions  .MAT, 
.INR,  and  .ABS  are  needed.  These  files  are  1)  an  ASCII  file  that  is  the  main  material 
file,  2)  an  ASCII  file  containing  wavelength  versus  reed  index  of  refraction  data,  and  3) 
an  ASCII  file  containing  wavelength  versus  absorption-coefficient  data,  respectively.  All  of 
the  file  names  are  stored  in  the  parameter  (PRM)  file.  In  default  operation,  the  INR  and 
ABS  files  will  have  the  same  name  as  the  MAT  file.  The  user  has  the  option  of  associating 
alternate  ABS  and  INR  files  with  a given  MAT  file  in  a particular  parameter  file,  even 
though  they  have  different  names.  If  this  is  done,  asterisks  follow  the  MAT  file  names  in  the 
“Materials:”  option  heading,  as  shown  in  figure  1.  Other  material  properties  are  included 
in  internal  models.  Some  of  these  allow  the  user  to  specify  the  values  of  the  parameters 
in  the  equations  defining  the  models,  and  most  can  be  replaced  by  user-supplied  ASCII 
files  having  the  appropriate  extensions.  The  internal  model  of  the  absorption  coefficient  of 
silicon,  which  is  discussed  next,  is  an  example. 

The  default  absorption-coefficient  data  in  the  SLABS  file  of  Version  2 of  PC-lD  are  not 
accurate  enough  for  high-accuracy  modeling  of  silicon  photodiodes  over  the  400  to  900  nm 
spectral  region  [3].  Therefore,  two  data  sets  of  higher  accuracy  were  calculated  and  stored 
in  files  named  SIL.WEAK.ABS  and  SILJPHIL.ABS.  The  first  data  set  was  calculated  by 
using  eq  (1)  of  Ref.  [3],  which  was  fitted  to  the  data  of  Ref.  [4]  over  the  spectral  range 
from  470  to  1180  nm.  The  second  data  set  was  calculated  by  interpolation  of  the  data 
of  Ref.  [5].  The  first  set  is  expected  to  be  more  accurate  at  longer  wavelengths,  and  the 
second  at  shorter  wavelengths. 

The  ratios  of  the  data  in  SIL_PHIL.ABS  to  that  in  SIL.WEAK.ABS  are  plotted  over  the 
400  to  900  nm  spectral  range  in  figure  2.  The  difference  can  be  characterized  as  12%  - 
4%  over  the  spectral  region  from  440  to  780  nm,  growing  much  larger  outside  that  spectral 
region.  The  value  of  12%  is  a convenient  average  offset  because  it  is  the  difference  between 
the  two  data  sets  at  633  nm.  At  that  wavelength,  the  value  calculated  from  eq  (1)  of  Ref. 
[3]  agrees  with  a recent  measurement  [6]  to  within  the  ±2%  uncertainty  associated  with 
the  measurement. 

One  material  property  not  modeled  by  PC-lD  is  the  quantum  yield  for  electron  hole  pair 
production  [7].  This  quantity  may  differ  significantly  from  unity  outside  the  400  to  900 
nm  spectral  region  [8],  but  no  high  accuracy  models  currently  exist  [7,9].  This  is  the  main 
reason  that  the  modeling  described  in  Part  II  of  this  paper  is  confined  to  the  400  to  900 
nm  spectral  region.  However,  other  problems,  such  as  larger  uncertainties  in  the  available 
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Figure  2.  Ratios  of  the  absorption-coefficient  data  interpolated  from  Ref.  [5]  to  those 
calculated  from  eq  (1)  of  Ref.  [3]  that  were  fitted  to  the  data  of  Ref.  [4]. 
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absorption-coefficient  data  at  shorter  wavelengths  and  uncertainties  in  the  fraction  of  the 
radiation  reflected  by  the  rear  surface  of  the  photodiode  at  the  longer  wavelengths,  also 
contribute  to  a significant  deterioration  in  accuracy  outside  that  spectral  region. 

2.1.2  Device  Properties 

The  “Doping:”  heading  of  the  opening  menu  of  PC-lD  allows  either  an  internal  doping 
model  or  an  external  file  containing  doping  concentrations  for  each  region  of  the  device 
being  modeled.  The  internal  model  consists  of  two  front  and  two  rear  dopant  distributions, 
as  well  as  a uniform  background  dopant.  The  distributions  can  be  chosen  from  uniform, 
complementary  error,  and  Gaussian  functions,  and  the  parameters  defining  the  maximum 
value,  its  location,  and  the  width  of  the  distributions  can  be  adjusted.  If  tabular  doping 
data  are  to  be  used,  they  are  read  from  a user-generated  ASCII  file  containing  a depth,  an 
n-type  dopant  concentration,  and  a p-type  dopant  concentration  on  each  line,  and  having 
.DOP  as  its  extension. 

The  shape  of  the  equilibrium  majority-carrier  concentration  near  the  front-surface  oxide- 
silicon  interface  is  very  important  to  high-accuracy  photodiode  modeling  [10].  It  is  neces- 
sary to  use  external  files  to  model  the  front  region  doping  in  p'^'n  type  photodiodes  to  force 
PC-lD  to  devote  enough  finite  elements  to  the  front  region  to  approximate  accurately  the 
majority- carrier  concentration  there. 

The  “Recombination:”  heading  of  the  opening  menu  of  PC-lD  allows  the  defect-related 
recombination  in  the  volume  and  at  the  surfaces  of  the  device  to  be  modeled.  (Auger 
recombination  is  considered  a material  property  and  is  covered  under  the  “Materials:” 
heading.)  Shockley- Read-Hall  (SRH)  recombination  through  a single  energy-level  trap 
state  is  used  as  the  model  for  volume  recombination.  Either  an  internal  model  or  a user- 
defined  table  of  depths  and  electron  and  hole  lifetimes  (equivalent  to  cross  sections  in  the 
SRH  model)  in  an  external  ASCII  file  with  the  extension  .TAU  can  be  used  to  model  the 
volume  recombination.  With  the  internal  model,  the  user  may  choose  a single  electron 
lifetime,  a single  hole  lifetime,  and  a single  trap  energy  relative  to  mid-gap  for  each  region. 
With  the  external  model,  a trap  level  at  mid-gap  is  used  with  the  lifetime  data  in  the 
external  ASCII  file. 

The  restriction  of  the  SRH  model  to  a single-energy  state  prevents  PC-lD  from  accurately 
fitting  the  measured  [11]  variation  of  quantum  efficiency  with  flux  (i.e.,  nonlinearity)  at 
950  nm  for  an  EG&G  UV444B  photodiode  as  illustrated  in  figure  3.  The  simulated  data 
in  figure  3 were  calculated  using  a mid-gap  state  with  equal  electron  and  hole  lifetimes. 
Even  though  the  simulated  and  experimental  data  agree  in  the  low  flux  limit  and  again  at 
a high  flux  level,  their  shapes  are  very  different.  Better  agreement  can  be  achieved  at  low 
flux  levels  at  the  cost  of  a worse  fit  at  the  higher  flux  levels  by  moving  the  state  away  from 
mid-gap,  and  by  using  different  electron  ajid  hole  lifetimes  as  illustrated  in  figure  4.  It  is 
expected  that  the  proper  distribution  of  SRH  trap  levels  over  the  band  gap  would  result 
in  an  accurate  simulation  of  the  nonlinearity  of  this  type  of  photodiode. 

It  should  be  possible  to  develop  an  external  lifetime  model  based  on  a parameterized  dis- 
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Figure  3.  Comparison  of  experimental  (filled  circles)  and  simulated  (filled  squares)  linearity 
measurements  on  an  EG<SiG  UV444B  photodiode  at  950  nm  using  a single  SRH  trap  level 
at  mid-gap  with  electron  and  hole  lifetimes  of  74.04  /xs. 
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Figure  4.  Comparison  of  experimental  (filled  circles)  and  simulated  (filled  squares)  linearity 
measurements  on  an  EG&G  UV444B  photodiode  at  950  nm  using  a single  SRH  trap  level 
0.1796  eV  above  mid-gap  with  electron  and  hole  (check  for  vice  versa)  lifetimes  of  900  ms 
and  700  ms,  respectively. 
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tribution  of  SRH  recombination  states,  and  to  adjust  the  parameters  to  fit  the  nonlinearity 
data  shown  in  figures  3 and  4.  However,  it  would  be  necessary  alternately  to  run  PC-lD 
to  calculate  the  carrier  concentrations  and  then  the  external  model  to  calculate  the  life- 
times appropriate  to  these  carrier  concentrations,  and  to  iterate  to  self-consistency.  This 
is  beyond  the  scope  of  the  work  reported  in  this  paper. 

The  failure  of  the  model  to  describe  the  nonlinearity  of  this  type  of  photodiode  casts  some 
doubt  on  its  ability  to  model  accurately  the  reverse-bias  self-calibration  experiment  :12-13\ 
This  point  is  addressed  in  more  detail  in  Part  II  of  this  paper. 

In  the  case  of  surface  recombination,  the  “Recombination:”  option  allows  the  choice  of 
one  of  two  internal  models,  a surface-recombination  velocity  model,  and  a saturation- 
current  density  model,  but  does  not  accommodate  a user-defined  data  file.  The  choice  of 
the  surface-recombination  velocity  model  allows  the  user  to  adjust  the  hole-  and  electron- 
recombination  velocities  as  well  as  the  energy-level  of  a single-energy  surface  state. 

The  “Surfaces:”  option  allows  the  user  to  choose  the  surface  charge  at  the  front  and  back 
surfaces  of  the  device.  For  modeling  the  oxide-bias  self- calibration  experiment  [12-I3y 
the  surface  charge  is  the  algebraic  sum  of  any  charge  trapped  in  the  oxide  as  a result  of 
thermal  oxidation  and  the  charge  stored  on  the  oxide  surface  by  the  voltage  applied  to  the 
transparent  electrode. 

The  “Circuit:”  option  of  the  main  menu  of  PC-lD  allows  the  user  to  define  a circuit  in 
which  the  device  is  to  operate.  For  the  case  of  a photodiode,  the  emitter  and  base  are 
connected,  the  collector  is  disconnected,  and  the  connections  to  the  device  are  made  at  the 
front  and  back  surface.  Keeping  the  collector  disconnected  forces  the  “Collector  Voltage:” 
option  to  zero,  as  shown  in  figure  1.  Setting  the  internal  resistance  to  some  non-zero 
value  allows  the  effect  of  series  resistance  to  be  simulated. 

3.  The  Photodiode  Modeling  Shell  Programs 

The  photodiode  modeling  shell  for  Version  2 of  PC-lD  is  designed  to  allow  simulation 
of  oxide-bias  and  reverse-bias,  self-calibration  experiments,  nonlinearity  measurements, 
and  internal  quantum  eflB.ciency  spectra,  while  accessing  the  highest  accuracy  and  preci- 
sion available  from  PC-lD.  It  consists  of  three  program  files.  The  first  program,  a short 
MS-DOS  Batch  program  named  RUN_PClD.BAT,  is  listed  in  figure  5.  The  other  two 
programs,  MAKE_PRM.EXE  and  READ  J^DF.EXE,  were  complied  with  Version  5.5  of 
Turbo  Pascal.  Listings  of  the  source  code  for  these  programs  are  given  in  Appendices  I 
and  II. 

3.1  Operation  of  RUN_PC1D.BAT 

The  first  thing  that  RUN_PC1D.BAT  does  when  run  is  to  test  whether  or  not  the  files 
MAKE-PRM.EXE,  PC-lD. EXE,  and  READJPDF.EXE  all  exist.  If  not,  it  issues  an  error 
message,  and  terminates  execution  at  :STOP.  If  all  three  files  do  exist,  RUN_PClD  deletes 
the  temporary  file  TEMP.DAT.  This  is  a precaution  in  case  a previous  nonstandard  ter- 
mination left  this  file  in  existence.  RUNJPClD  then  enters  the  loop  between  :LOOP  and 


10 


IF  NOT  EXIST  MAKE_PRM.EXE  GOTO  ERROR  1 
IF  NOT  EXIST  PC-1D.EXE  GOTO  ERROR2 
IF  NOT  EXIST  READ_PDF.EXE  GOTO  ERRORS 
DEL  TEMP.DAT 
: LOOP 

MAKE_PRM.EXE 

IF  NOT  EXIST  TEMP.DAT  GOTO  STOP 

PC-1D  TEMP  TEMP 

READ__PDF  . EXE 

GOTO  LOOP 
: ERROR1 

REM  MAKE__PRM.EXE  not  found. 

GOTO  :MESSAGE 
: ERROR2 

REM  PC-1D.EXE  not  found. 

GOTO  MESSAGE 
: ERRORS 

REM  READ_PDF.EXE  not  found. 

GOTO  MESSAGE 
:MESSAGE 

TYPE  RUNPC1D.DOC 

GOTO  STOP 
:STOP 


Figure  5.  Listing  for  the  MSDOS  batch  file  RUNJPC1D.BAT  that  supervises  the  execution 
of  the  programs  MAKE_PRM.EXE,  PC-lD.EXE,  and  READ_PDF.EXE  for  high-accuracy 
photodiode  modeling. 
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GO  TO  LOOP.  Once  in  this  loop,  it  transfers  control  to  MAKE_PRM.EXE. 

The  first  thing  that  MAKE_PRM.EXE  does  is  to  test  the  existence  of  TEMP.DAT.  If 
TEMP.DAT  does  not  exist,  then  MAKE_PRM.EXE  prompts  the  user  to  define  the  pho- 
todiode experiment  to  be  simulated.  Examples  of  the  questions  and  typical  answers  are 
shown  in  Appendix  A.  The  questions  and  answers  are  straightforward  and  require  no 
explanation.  If  the  file  option  rather  than  the  keyboard  option  is  chosen  to  define  the 
independent  variables  for  the  experiment  to  be  simulated,  then  the  program  asks  for  the 
file  name  instead  of  the  start,  stop,  and  step  values  for  the  independent  variable.  Once 
MAKE_PRM.EXE  has  obtained  all  of  the  necessary  information,  it  writes  that  information 
into  the  file  TEMP.DAT,  creates  the  temporary  parameter  file  TEMP.PRM,  and  passes 
control  back  to  RUN_PC1D.BAT. 

RUN_PC1D.BAT  now  tests  the  existence  of  TEMP.DAT.  If  it  does  not  exist,  RUN_PClD 
terminates  execution.  If  TEMP.DAT  does  exist,  then  RUN_PC1D  passes  control  to  PC- 
ID.EXE.  This  program  reads  TEMP.PRM  (the  first  TEMP  on  the  command  line  following 
PC-lD),  carries  out  the  simulation  defined  therein,  writes  a table  having  a depth,  an 
electron  current,  a hole  current,  and  a total  current  on  each  line  in  the  file  TEMP.PDF 
(the  second  TEMP  on  the  command  line),  and  returns  control  to  RUN_PClD,BAT.  Note 
that  PDF  files  are  the  standard  ASCII  format  HP  plotter  output  files  generated  by  PC-lD. 

RUN_PC1D.BAT  now  passes  control  to  READ_PDF.EXE.  This  program  reads  the  total 
photocurrent  at  the  depth  specified  in  the  definition  of  the  experiment  to  be  simulated 
(see  Appendix  A)  from  file  TEMP.PDF,  and  appends  it  to  the  output  file  specified  in 
response  to  the  prompt  from  MAKE_PRM.EXE.  READ_PDF.EXE  then  returns  control 
to  RUxN_PClD.BAT. 

On  subsequent  passes  through  the  loop,  MAKE_PRM.EXE  finds  that  TEMP.PRM  does 
exist,  reads  the  necessary  information  from  this  file,  and  writes  updated  TEMP.DAT  and 
TEMP.PRM  files.  The  information  about  when  to  exit  the  loop  is  in  TEMP.DAT,  and 
when  MAKE_PRM.EXE  determines  that  this  time  has  come,  it  deletes  TEMP.DAT,  trig- 
gering RUN_PC1D.BAT  to  terminate  execution. 

3.2  Features  of  MAKE_PRM  and  READ_PDF 

A few  points  about  MAKE_PRM.EXE  and  READ_PDF.EXE  are  in  order.  First,  these  two 
programs  are  designed  to  compensate  for  traits  of  PC-lD  that  make  it  less  than  ideal  for 
high-accuracy  modeling  of  photodiodes.  The  first  of  these  is  that  PC-lD. EXE  calculates 
the  internal  quantum  efficiency  from  the  sum  of  the  photocurrent  and  the  dark  current 
occurring  under  reverse  or  forward  bias.  This  may  not  be  a problem  in  solar  cell  modeling, 
but  it  is  not  correct  for  photodiode  modeling.  As  a resiilt,  MAKE_PRM.EXE  calculates 
and  stores  the  appropriate  dark  currents  for  all  simulations  that  have  a nonzero  bias 
voltage  Vfcb,  and  READ_PRM.EXE  reads  them  and  subtracts  them  from  the  appropriate 
total  currents  when  calculating  the  quantum  efficiency.  This  is  simulated  as  a shutter- 
closed/shutter-open  measurement  sequence  by  MAKE_PRM.EXE. 

Even  though  the  total  current  is  a constant  throughout  a photodiode,  it  must  be  calculated 
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at  each  finite  element  in  the  photodiode.  The  numerical  precision  with  which  it  can  be 
calculated  depends  upon  the  majority-carrier  concentration  in  the  element  where  it  is  being 
calculated  due  to  nature  of  the  numerical  algorithm,  which  involves  the  difference  of  two 
large  numbers  that  scale  with  doping  concentration.  This  is  illustrated  in  figure  6,  which 
shows  the  variation  in  the  total  current  computed  by  PC-lD  as  a function  of  depth  in  the 
front  region  of  a p^n  type  photodiode  irradiated  by  1 /iW/cm^  of  440  nm  radiation.  This 
result  was  obtained  by  running  PC-lD  through  its  standard  user  interface  and  could  not 
be  obtained  by  running  PC-lD  with  the  shell  just  described,  since  carrier  concentration 
as  a function  of  depth  is  not  an  output  supported  by  the  shell.  The  decrease  in  the  noise 
in  the  data  with  increasing  distance  from  the  front  surface  of  the  photodiode  is  evident  in 
the  figure. 

PC-lD. EXE  calculates  the  quantum  efficiency  from  the  total  current  at  either  the  front  or 
the  rear  of  the  photodiode  at  the  discretion  of  the  user.  Since  these  are  usually  the  regions 
where  the  majority-carrier  concentration  is  greatest,  PC-lD  often  calculates  the  quantum 
efficiency  from  values  of  the  total  photocurrent  that  are  not  the  most  precise  available  to 
it.  This  loss  of  precision,  which  occurs  even  in  the  absence  of  dark  current,  is  the  reason 
that  MAKE_PRM.EXE  asks  the  user  to  choose  the  depth  at  which  the  total  photocurrent 
will  be  read  from  the  TEMP.PDF  file.  This  depth  should  be  somewhere  in  the  interior  of 
the  photodiode  where  the  doping  concentrations  are  low  and  the  effect  of  surface  fields  on 
the  carrier  concentrations  is  small  or  negligible. 

Very  high  precision  is  available  under  these  conditions.  The  internal  quantum  efficiency 
calculated  by  RUN_PClD  for  the  same  type  of  photodiode  used  to  generate  the  data  in 
figure  6,  but  with  aU  loss  (recombination)  mechanisms  set  to  zero,  is  6 parts  per  million 
(ppm)  above  unity  ±1  ppm  for  1 /iW/cm^  irradiation  anywhere  within  400  to  900  nm 
spectral  region.  This  result  was  obtained  by  assuming  the  somewhat  large  thickness  of  400 
mm  for  this  photodiode  to  minimize  the  fraction  of  the  incident  radiation  penetrating  to 
the  rear  of  the  photodiode.  Since  PC-lD  does  not  permit  a reflectance  of  100%  at  the  rear 
surface  of  the  diode,  a more  typical  300-mm  thick  photodiode  would  show  a loss  of  a few 
ppm  at  the  longest  wavelengths. 

Since  PC-lD  does  not  model  the  oxide  passivation  on  the  front  surface  of  the  photodiode, 
it  is  necessary  to  vary  the  oxide  trapped  charge  N,,  to  simulate  the  effect  of  applying  a 
voltage  across  the  oxide  in  the  silicon  self-calibration  experiment.  First,  the  value  of  N,, 
must  be  chosen  to  simulate  the  oxide-fixed  chaxge.  To  this  value  must  be  added  the  charge 
appropriate  to  the  application  of  oxide  bias  according  to  the  formula  for  the  charge  density 
stored  on  a capacitor  as  a function  of  thickness,  dielectric  constant,  and  applied  voltage. 
It  is  this  sum  that  PC-lD  must  use  as  N,,  when  simulating  the  effect  of  oxide  bias  on 
a photodiode.  The  bookkeeping  for  the  values  of  N,,  chosen  to  simulate  the  oxide-fixed 
charge  and  the  values  of  N,,  needed  to  simulate  the  combined  effect  of  the  fixed  charge 
and  the  applied  oxide  bias  are  handled  by  MAKE-PRM.EXE. 
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DISTANCE  FROM  FRONT  (|im) 

Figure  6.  Variation  in  the  total  photocurrent  calculated  by  PC-lD  as  a function  of  position 
in  the  front  region  of  a Hamamatsu  1337  type  photodiode  due  to  the  effect  of  the  high 
front  region  doping  concentration. 
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4.  Conclusion 


The  semiconductor  device  modeling  program  PC-lD  and  some  programs  to  enhance  its  use 
in  the  high-accuracy  modeling  of  silicon  photodiodes  have  been  described.  The  features  of 
PC-lD  that  suit  it  for  high-accuracy  photodiode  modeling  have  been  described,  and  the 
features  of  the  support  programs  that  compensate  for  features  of  the  PC-lD  user  interface 
that  are  less  than  ideal  for  this  application  have  also  been  described.  Examples  of  the  use 
of  these  programs  in  photodiode  modeling  in  different  types  of  high-accuracy  applications 
are  given  in  Parts  II  and  III  of  this  paper. 
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APPENDIX  I.  A listing  of  MAKE_PRM.PAS  that  compiles  under  Turbo  Pascal  Version  5.5 
PROGRAAl  mak:e_temporary_prm_files_for_running_pc  1 d; 

{RUN_PC1D.BAT  runs  this  program  when  this  program  is  named  MAKE_PRM.EXE. 

This  program  writes  2 temporary  files,  TEMP.PRM  and  TEMP.DAT. 

RUN_PC1D.BAT  then  runs  PC-ID  with  TEMP.PRM  as  input  and  TEMP.PDF  as  output. 
RUN_PC1D.BAT  then  runs  READ_PDF.EXE. 

READ_PDF.EXE  uses  TEMP.PDF  and  TEMP.DAT  to  write  QE  data  to  a file. 
RUN_PC1D.BAT  runs  MAKE_PRM.EXE  again. 

MAKE_PRM.EXE  uses  TEMP.DAT  to  decide  to  either 

1)  increment  some  parameters  in  TEMP.DAT  and  TEMP.PRM,  or 

2)  delete  the  3 temporary  files.} 

{$N+} 

USES  CRT,  DOS; 

TYPE  file_of_values_type  = FILE  OF  DOUBLE; 

TYPE  pjndex_type  = 

(I_parameters,  I_maxJterations, 

I_mode,  I_graphs,  I_size_factor,  I_error Jimit, 

I_clamp,  I_step_size,  I_final_step, 

I_solutions_step,  I_eq_renode,  I__ss_renode,  I_tr_renode, 

I_graph_l,  I_graph_2,  I_graph_3,  I__graph_4, 

I_area,  I_area_value,  I_facet_depth,  I_facet_angle, 

I_temperature,  I_temp_units, 

I_front_surface,  I__back_surface, 

I_generation,  I_front_texture,  I_back_texture, 

I_int__refl_value,  I_front_refl, 

I_front_refl_value,  I_back_refl_value, 

I_front_nl,  I_front_dl, 

I_front_n2,  I_front__d2, 

I_front_n3,  I_front_d3, 

I_front_barrier,  I_back_barrier, 

I_front_Nss,  I_back_Nss, 

I_base_location,  I_emitter_location,  I_collectorJocation, 

I_shunt_locl,  I_shunt_loc2, 

I_band_file,  I_boundary__dx,  I_recombine, 

I_collector,  I_base,  I_shunt, 

I_Vcc_shape,  I_Vcc_start,  I_Vcc_end, 

I_Vbb_shape,  I_Vbb_start,  I_Vbb_end, 

I_light_shape,  IJight_start,  I_light_end, 

I_spectrum,  I_lambda_start,  I_lambda_end, 

I_Rcc,  I_Rbb,  I_Rc,  I_Rb, 
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I_Gb,  I_Gc,  I_Gsh, 

I_Db,  I_Dc,  I_Nb,  I_Nc); 

q_index_type  = 

(I_w,  I_material,  I_material_mod, 

I_doping,  I_recomb_front,  I_recomb_back;, 

I_taun,  I_taup,  I_Et, 

I_Sn_front,  I_Sp_front,  I_Et_front, 

I_Sn_back,  I_Sp_back,  I_Et_baclc, 

I_JO_front,  I_JO_front_n,  I_JO_back,  I_JO_back_n, 

I_Nbk,  I_Nbk_type, 

I_Nfl_type,  I_Nfl_shape,  I_Nfl_depth,  I_Nfl_peak, 
I_Nf2,  I_Nf2_type,  I__Nf2_shape,  I_Nf2_depth,  I_Nf2_peak, 
I__Nrl,  I_Nrl_type,  I__Nrl_shape,  I_Nrl_depth,  I_Nrl_peak, 
I_Nr2,  I_Nr2_type,  I_Nr2_shape,  I_Nr2_depth,  I_Nr2_peak, 
I_Eg,  I_affinity,  I_NcNv,  I_niO,  I_eps, 

I_Cnt,  I_Cpt,  I_Cn,  I_Cp,  I_Bn,  I_Bp, 

I_BGN_Phid,  !_BGN_Nd,  I_BGN_Phia,  I_BGN_Na, 
I_n_rile, 

I_n, 

I_mu_max_n, 

I_mu_min_nn,  I_Nref_nn,  I_alpha_nn, 

I_mu_min_np,  I__Nref_np,  I_alpha_np, 

I_mu_nbl,  I_mu_nb2,  I_mu_nb3,  i_mu_nb4, 

I_mu_max_p, 

I_mu_min_pp,  I_Nref_pp,  I_alpha_pp, 

I_mu_min_pn,  I_Nref_pn»  I_alpha_pn, 

I_mu_pbl,  I_mu_pb2,  I_mu_pb3,  i_mu_pb4, 

I_n_sat,  I_p_sat, 

I_absorb_file, 

I_Egdl,  I_Egd2,  I_Eg_l,  I_Eg2,  I_Epl,  I_Ep2, 

I_Asb_Adl,  I_Abs_Ad2, 

I_Abs_Al,  I_Abs_A2,  I_Abs_A3, 1_Abs_A4, 

I_Abs_B,  I_Abs_Tg); 

variant_type  = (str_16,  byt_16,  rea_08,  int__04); 

s_type  = STRING[15]; 
b_type  = ARRAY[0..15]  OF  BYTE; 
r_type  = ARRAY[0..1]  OF  DOUBLE; 
i_type  = ARRAY[0..3]  OF  LONGINT; 

p_type  = RECORD 

CASE  variant  : variant_type  OF 
str_16  : (s  : ARRAY[pJndex_type]  OF  s_type); 
byt_16  : (b  : ARRAY[p_index_type]  OF  b_type); 
rea_08  : (r  : ARRAY[p_index__type]  OF  r_type); 
int_04  : (i  : ARRAY[pJndex_type]  OF  i_type); 
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END; 


q_type  = RECORD 

CASE  variant  : variant_type  OF 

str_16  : (s  : ARRAY[1..3]  OF  ARRAY [q_index_type]  OF  s_type); 
byt_16  : (b  : ARRAY[1..3]  OF  ARRAY[q_index_type]  OF  b_type); 
rea_08  : (r  : ARRAY[1..3]  OF  ARRAY[c|_index_type]  OF  r_type); 
int__04  : (i  : ARRAY[1..3]  OF  ARRAY[q_index_type]  OF  i_type); 
END; 

TYPE  prm_data_type 
= OBJECT 

model_pnn_rile_path  : PATHSTR; 
p : p.type; 
q : q_type; 

PROCEDURE  read_j>nn_file;  {page  > 

PROCEDURE  write_prm_file;  {page  } 

END; 

TYPE  exp_data_type 

= OBJECT(prm_data_type) 
key  : CHAJl; 
shutter  : STRING[6]; 

lambda,  energy,  phi,  Vo,  Vr,  depth,  oxide_t,  value  : DOUBLE; 
index,  last  : LONGINT; 

QE_file_path  : PATHSTR; 

FUNCTION  file_exists(  path  : PATHSTR)  : BOOLEAN;  {page  } 


PROCEDURE  initialize_or_update;  {page  } 

PROCEDURE  initialize;  {page  } 

PROCEDURE  get_parameter_to_scan;  {page  } 

PROCEDURE  get_values_from_file;  {page  } 

{OR} 

PROCEDURE  get_start_stop_and_step;  {page  } 

PROCEDURE  setj)arameter_to_value;  {page  } 

PROCEDURE  open_empty_QE_file;  {page  } 

PROCEDURE  write_temp_data;  {page  } 

PROCEDURE  update;  {page  } 

PROCEDURE  get_current_value; 

{ set_parameter_to_value; } 

{ write_temp_data; } 

PROCEDURE  prepare_pq;  {page  } 

PROCEDURE  delete_temporary_files;  {page  } 

END; 


VAR  exp_data  : exp_data_type; 

FUNCTION  exp_data_type.file_exists(  path  : PATHSTR)  : BOOLEAN; 
VAR  f : TEXT; 
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BEGIN 

{$!-} 

ASSIGN(f,  path); 

RESET(f); 

CLOSE(0; 

{$!+} 

IF  (lORESULT  = 0) 

THEN  file_exists  :=  TRUE 
ELSE  file_exists  :=  FALSE; 

END; 

PROCEDURE  exp_data_type.initialize_or_update; 

VAR  f : TEXT; 

BEGIN 

IF  NOT  file_exists(  ’temp.dat’ ) 

THEN  BEGIN 
CLRSCR; 

IF  file_exists(  ’JO.dat’  ) 

THEN  BEGIN 

ASSIGN(f,  ’JO.dat’  ); 

ERASE(f); 

END; 

WRITELN( 

’Parameter  file  (including  PRM  extension)  defining  photodiode:  ’); 
READLN(modelj)rm_file_path); 

IF  NOT  file_exists(  model_prm_file_path  ) 

THEN  BEGIN 

WRlTELN(’File:  ’,  modelj)rm_file_path,  ’ not  found.  Halting.’); 
HALT; 

END; 

WRITELN( ’Complete  path  to  file  in  which  to  store  QE  results:  ’); 

READLN(QE_file_path); 

initialize; 

END 

ELSE  BEGIN 
update; 

END; 

END; 

PROCEDURE  exp_data_type.initialize; 

VAR  s : STRING[1]; 

BEGIN 

get_parameter_to_scan; 
s :=  ’ ’; 

WHILE  NOT  ( UPCASE(s[l])  IN  [’F,  ’K’]  ) 

DO  BEGIN 
WRITE( 

’Start,  step,  and  stop  values  from  (K)eyboard  or  all  values  from  (F)ile?  ’); 
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READLN(s); 
s :=  UPCASE(s[l]); 

IF  NOT  ( UPCASE(s[l])  IN  [’K’,  ’F]  ) 

THEN  WRrTELN(’Your  response  is  unacceptable.  Try  it  again.’); 

END; 

IF  s = ’F’ 

THEN  get_values_from_file 
ELSE  get_start_stop_and__step; 

IF  NOT  (key  IN  [T’,  ’W’]) 

THEN  BEGIN 

WRITE( ’Wavelength  (nm):  ’); 

READLN(lambda); 

END; 

IF  key  <>  ’L’ 

THEN  BEGIN 

WRrTE( ’Incident  irradiance  (W  per  sq.  cm):  ’); 

READLN(phi); 

END; 

IF  key  <>  ’O’ 

THEN  BEGIN 

WRITE(’Oxide  bias  voltage  (V):  ’); 

READLN(Vo); 

END; 

IF  key  < > ’R’ 

THEN  BEGIN 

WRITE( ’Reverse  bias  voltage  (V):  ’); 

READLN(Vr); 

END; 

WRITE( ’Non-zero  thickness  of  front  surface  oxide  (nm):  ’);  READLN(oxide_t); 
WRITEfDepth  at  which  to  get  total  photocurrent  (um):  ’);  READLN(depth); 
IF  Vr  = 0 

THEN  shutter  :=  ’open’ 

ELSE  shutter  :=  ’closed’; 

open__empty_QE_file; 

write_temp_data; 

END; 


PROCEDURE  exp_data_type.get_parameter_to_scan; 

TYPE  letter_set_type  = SET  OF  ’A’  ..  ’Z’; 

CONST  answer_set  : letter_set_.type  = [’O’,  ’R’,  ’W’,  ’P’,  ’L’]; 
VAR  s : STRING[1]; 

BEGIN 
s :=  ’ ’; 

WHILE  NOT  ( UPCASE(s[l])  IN  answer_set ) 

DO  BEGIN 
{CLRSCR;} 

WRITELN(Type  of  experiment  to  simulate?  ’); 
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WRITELN(’(0)xide  bias  (Volts)  scan?  ’); 

WRITELN(’(R)everse  bias  (Volts)  scan?  ’); 

WRITELN(’(W)avelength  (nm)  scan?  ’); 

WRITELN(’(P)hoton  energy  (eV)  scan?  ’); 

WRITELN(’(L)inearity  (W/cm^2)  scan?  ’); 

WRITE(’Enter  letter  key  here:  ’);  READLN(s); 

IF  NOT  ( UPCASE(s[l])  IN  answer_set  ) 

THEN  \WITELN(’Your  response  is  unacceptable.  Try  it  again.’); 

END; 

key  :=  UPCASE(s[l]); 

END; 

PROCEDURE  exp_data_type.get_values_from_file; 

VAR  start,  stop,  step  : DOUBLE; 
name  : STRING[80]; 
f : TEXT; 

g : file_of_values_type; 

BEGIN 

WRrrE(’Name  of  file  with  all  values  of  scanned  parameter:  ’);  READLN(name); 
ASSIGN(f,  name); 

RESET(f); 

ASSIGN(g,  ’temp.VAL’); 

REWRrTE(g); 

READLN(f,  start); 

WRITE(g,  start); 
last  :=  0; 

WHILE  NOT  EOF(0 
DO  BEGIN 
READLN(f,  stop); 

WRrTE(g,  stop); 

INC(last); 

END; 

WRrrE(g,  stop); 

INC(last); 

CLOSE(f); 

CLOSE(g); 

index  :=  0; 

value  :=  start; 

set_parameter_to_value; 

END; 

PROCEDURE  exp_data_type.get_start_stop_and_step; 

VAR  start,  stop,  step  : DOUBLE; 

g : file_of_values_type; 

BEGIN 

WRi  rE(’First  value  for  scanned  parameter:  ’);  READLN(start); 

WRITEfStop  value  for  scanned  parameter:  ’);  READLN(stop); 

IF  key  = ’L’ 
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THEN  BEGIN 

WRrrE(’Factor  by  which  to  increase  irradiance  on  each  iteration:  ’); 
READLN(step); 

IF  ( (start>stop)  AND  (step>  = l)  ) 

OR  ( (start < stop)  AND  (step<“l)  ) 

THEN  BEGIN 

WRITELN( These  values  would  produce  an  endless  loop.  Halting.’); 
HALT; 

END; 

END 

ELSE  BEGIN 

WRITE(’Step  size  for  scanned  parameter:  ’); 

READLJ^I(step); 

IF  ( (start > stop)  AND  (step>=0)  ) 

OR  ( (start<stop)  AND  (step<-0)  ) 

THEN  BEGIN 

WRrTELN( These  values  would  produce  an  endless  loop.  Halting.’); 
HALT; 

END; 

END; 

ASSIGN(g,  ’temp.VAL’); 

REWRITE(g); 
last  :=  0; 

IF  key  - L’ 

THEN  BEGIN 

value  start/step; 

IF  start  < stop 

THEN  WHILE  value  < - stop 
DO  BEGIN 
value  :=  value* step; 

WRITE(g,  value); 

INC(last); 

END 

ELSE  WHILE  value  > = stop 
DO  BEGIN 
value  :=  value* step; 

WRITE(g,  value); 

INC(last); 

END; 

END 

ELSE  BEGIN 

value  :=  start-step; 

IF  start  < stop 

THEN  WHILE  value  <=  stop 
DO  BEGIN 
value  valued  step; 

WRITE(g,  value); 

INC(last); 
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END 

ELSE  WHILE  value  >=  stop 
DO  BEGIN 
value  :=  value  + step; 

WRrTE(g,  value); 

INC(last); 

END; 

END; 

CLOSE(g); 

DEC(last); 

index  ;=  0; 

value  :=  start; 

set_parameter_to_value; 

END; 

PROCEDURE  exp_data_type.set_parameter_to_value; 

BEGIN 

CASE  key  OF 

’O’  : Vo  :=  value; 

’R’  : Vr  :=  value; 

’W’  : lambda  :=  value; 

’P’  : BEGIN 

energy  :=  value; 
lambda  :=  1239.85 /energy; 

END; 

’L’  : phi  :=  value; 

END;  {CASE} 

END; 

PROCEDURE  exp_data_type.open_empty_QE_file; 
VAR  f ; TEXT; 

BEGIN 

ASSIGN(f,  QE_file_path); 

REWRITE(f); 

CLOSE(f); 

END; 

PROCEDURE  exp_data_type.write_temp_data; 

VAR  f : TEXT; 

BEGIN 

ASSIGN(f,  ’temp.dat’); 

REWRITE(f); 

WRITELN(f,  model_prm_file_path); 

IF  key  <>  ’P’ 

THEN  WRITELN(f,  key,  ’ ’,  lambda,  ’ ’,  phi) 

ELSE  WRITELN(f,  key,  ’ ’,  energy,  ’ ’,  phi); 
WRITELN(f,  Vo,  ’ ’,  Vr,  ’ ’,  shutter); 

WRITELN(f,  depth,  ’ ’,  oxide^t); 
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WRITELN(f,  QE_file_path); 

WRrrELN(f,  index,  ’ last); 

WRITE LN( index,  last); 

CLOSE(f); 

END; 

PROCEDURE  exp_data_type.update; 

VAR  dir  : DIRSTR; 
name  : NAMESTR; 
ext  : EXTSTR; 
f : TEXT; 

BEGIN 

ASSIGN(f,  ’temp.dat’); 

RESET(0; 

READLN(f,  model_prm_file_path); 

IF  key  <>  T’ 

THEN  READLN(f,  key,  lambda,  phi) 

ELSE  BEGIN 

READLN(f,  key,  energy,  phi); 
lambda  :=  1239.85/energy; 

END; 

READLN(f,  Vo,  Vr,  shutter); 

WHILE  shutter[l]  = ’ ’ 

DO  DELETE(shutter,  1,  1); 

READLN(f,  depth,  oxide_t); 

READLN(f,  QE_file_path); 

READLN(f,  index,  last); 

CLOSE(f); 

IF  shutter  = ’open’ 

THEN  BEGIN 
INC(  index); 

IF  key  = ’R’ 

THEN  shutter  :=  ’closed’; 

END 

ELSE  shutter  :=  ’open’; 
get_current_value; 
set_parameter_to_value; 
write_temp_data; 

END; 

PROCEDURE  exp_data_type.get_current_value; 
VAR  g : file_of_values_type; 

BEGIN 

ASSIGN(g,  ’temp.val’); 

RESET(g); 

SEEK(g,  index); 

READ(g,  value); 

CLOSE(g); 
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END; 


PROCEDURE  exp_data_type.prepare__pq; 

CONST  eps_Si02  = 3.9; 
eps_Si  = 11.9; 

epsO  = 8.8541878e-14;  {C/(V  cm)} 
cm_per_nm  = le-7; 
q_per_C  = l/1.60218e-19; 

{and  from  W.  R.  Thurber  and  J.  R.  Lowney,  "Electrical  Transport 
Properties  of  Silicon"  Chapter  14  in  the  VLSI  Handbook,  ed.  N. 
G.  Einspruch  (Academic  Press,  Orlando,  1985),  page  188.} 

Ec  = 30.5e4; 
c = 0.657; 

VAR  NssO,  Nss,  RbO,  Rb,  EnO,  En,  bot,  top, 
sheet_R_ratio,  concentration_ratio  : DOUBLE; 

BEGIN 

IF  shutter  = ’closed’ 

THEN  phi  :=  0; 
p.s[I__mode]  :=  ’Steady  State’; 
p.s[I_graphs]  :=  ’Specified’; 
p.i[I _jraph_l][0]  :=  5; 
p.s[I_front_surface]  :=  ’Charged’; 
p.s[I_spectrum]  :=  ’Fixed’; 
p.s[I_light_shape]  :=  ’Fbted’; 
p.s[I_Vbb_shape]  :=  ’Fixed’; 
p.r[IJambda_start][0]  lambda; 
p.r[IJight_start][0]  :=  phi; 
p.r[I_Vbb_start][0]  :=  Vr; 

NssO  :=  p.r[I__front_Nss][0]; 

Nss  :=  NssO  + q_per_C*Vo*eps_SiO2*eps0/(oxide_t*cm_per_nm); 
p.r[I_front_Nss][0]  :=  Nss; 

RbO  :=  p.r[I_Rb][0]; 

IF  ( (Nss0=0)  OR  (Rb0=0)  ) 

THEN  BEGIN 
Rb  :=  RbO; 

END 

ELSE  BEGIN 
IF  Nss  = 0 

THEN  Nss  :=  Nss  + le7; 
concentration_ratio  :=  NssO/Nss; 

EnO  :=  ABS(NssO/(2*eps_Si*epsO)); 

En  :=  ABS(Nss/(2*eps_Si*eps0)); 
top  :=  1 -I-  EXP(c*LN(En/Ec)); 
bot  :=  1 + EXP(c*LN(EnO/Ec)); 
sheet_R_ratio  :=  top/bot; 

Rb  :=  RbO*concentration_ratio*sheet_R_ratio; 

WRlTELN(’conc.  ratio  = ’,  concentration_ratio); 
WRITELN(’sheet  R ratio  = ’,  sheet_R_ratio); 
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END; 

p.r[I__Rb][0]  :=  Rb; 

END; 

PROCEDURE  exp_data_type.delete__temporary_riles; 
VAR  f : FILE; 
dir  : DIRSTR; 
name  : NAMESTR; 
ext  : EXTSTR; 
temp_rile_path  : PATHSTR; 

BEGIN 

ASSIGN(f,  ’temp.dat’); 

ERASE(f); 

FSPLrr( model jprm_file__path,  dir,  name,  ext); 
temp_file_path  :=  dir  + ’temp.val’; 

ASSIGN(f,  temp__file_path); 

ERASE(f); 

temp_file_j)ath  :=  dir  + ’temp.prm’; 

ASSIGN(f,  temp__rile__path); 

ERASE(f); 

temp_file_j)ath  :=  dir  + ’temp.pdf; 

ASSIGN(f,  temp_file_path); 

ERASE(f); 

END; 

PROCEDURE  prm_data_typ€.readj)rm_rile; 

VAR  f : FILE; 

count,  result  : WORD; 

BEGIN 

ASSIGN(f,  model_prm__fiIe_path); 

RESET(f,  1); 

count  :=  16*SUCC(ORD(LNc)); 

BLOCKREAD(f,  exp_data,p.s,  count,  result); 
count  :=  16*3*SUCC(ORD(LAbs__Tg)); 
BLOCKREAD(f,  exp^data.q.s,  count,  result); 
CLOSE(f); 

END; 

PROCEDURE  prm_data_type.write_j)rm_file; 

VAR  f : HLE; 
dir  : DIRSTR; 
name  : NAMESTR; 
ext  : EXTSTR; 

temp_j)rm_file_path  : PATHSTR; 
count,  result  : WORD; 

BEGIN 

FSPLIT(model__prm__rile_j)ath,  dir,  name,  ext); 
temp_prm__file_path  :=  dir  + ’temp.prm’; 
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ASSIGN(f,  temp_prm_rile_j)ath); 

REWRITE(f,  1); 

count  :=  16*SUCC(ORD(I_Nc)); 
BLOCKWRrrE(f,  exp_data.p.s,  count,  result); 
count  :=  16*3*SUCC(ORD(I_Abs_Tg)); 
BLOCKWRlTH(f,  exp_data.q.s,  count,  result); 
CLOSE(f); 

END; 


PROCEDURE  test_exp_data; 

BEGIN 

WITH  exp_data 
DO  BEGIN 

WRITELN(model_prm_file_path); 

WRITELN(p.s[I_mode]); 

WRITELN(p.s[I_graphs]); 

WRrrELN( ’Graph  number  = ’,  p.i[I_SJ'aph_l][0]); 
WRITE  LN(p.s[I_spectrum]); 
WRITELN(p.s[I_light_shape]); 

WRrTELN( ’lambda  = ’,  p.r[l_lambda_start][0]); 
WRITELNC’phi  = ’,  p.r[IJight_start][0]); 
WRITELN(’Vo  = ’,  Vo); 

WRITELN(’Nss  = ’,  p.r[I_front_Nss][0]); 
WRITELN(’Vr  = ’,  p.r[I_Vbb_start][0]); 

END; 

END; 


BEGIN  {main 


WITH  exp_data 
DO  BEGIN 

WRITELN( ’MAKE_PRM.EXE  is  initializing  or  updating  as  appropriate.’); 
initialize_or__update; 

WRITELN(’The  shutter  is  now  ’,  shutter); 

IF  ( index  > = last ) 

THEN  BEGIN 

WRrTELN( ’MAKE_PRM.EXE  is  deleting  the  temporary  files.’); 
delete_temporary_files; 

END 

ELSE  BEGIN 

WRITELN(’MAKE_PRM.EXE  is  reading  ’,  model_prm_file_path); 
read_prm_file; 

{test_exp_data;} 

WRITELN(’MAKE_PRM.EXE  is  preparing  the  temporary  parameter  file.’); 
prepare_pq; 

{test_exp_data;} 

WRITELN(’MAKE_PRM.EXE  is  writing  TEMP.PRM.’); 
write_prm_file; 
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WRITELN; 

END; 

END; 

END. 
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APPENDIX  II.AListing  of  READ_PDF.PAS  that  compiles  under  Turbo  Pascal  Version  5.5. 
PROGRAlM  read_total_photocurrent_from_pdf_file; 

{$N  + } 

USES  CRT,  DOS; 

TYPE  pdf_data_type 
= OBJECT 

modelprm  filepath,  QE  filepath  : PATHSTR; 

key  : CHAU; 

shutter  : STRING[7]; 

lambda,  energy,  phi,  Vo,  Vr, 

oxide_t,  depth,  JO,  Jt,  QE  : DOUBLE; 

z : ARRAY[0..500]  OF  STRING[30]; 

PROCEDURE  initialize; 

PROCEDURE  read_pdf_file; 

PROCEDURE  write_JO_to_file; 

PROCEDURE  read_JO_from_file; 

PROCEDURE  append_QE_file; 

END; 

VAR  pdf_data  : pdf_data_type; 

PROCEDURE  pdf_data_type.initialize; 

VAR  f : TEXT; 

BEGIN 

ASSIGN(f,  ’temp.dat’); 

RESET(f); 

READLN(f,  model_prm_file_path); 

READLN(f,  key,  lambda,  phi); 

IF  key  = ’P’ 

THEN  BEGIN 

energy  :=  lambda; 
lambda  :=  1239.85/energy; 

END; 

READLN(f,  Vo,  Vr,  shutter); 

WHILE  shutter[l]  = ’ ’ 

DO  DELETE(shutter,  1,  1); 

READLN(f,  depth,  oxide_t); 

READLN(f,  QE_file_path); 

CLOSE(f); 

. END; 

PROCEDURE  pdf_data_type.write__JO_to_file; 

VAR  f : TEXT; 

BEGIN 
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{$!-} 

ASSIGN(f,  ’JO.dat’); 

RESET(0; 

CLOSE(0; 

{$!+} 

ASSIGN(f,  ’JO.dat’); 

IF  ( lORESULT  = 0 ) 

THEN  APPEND(f) 

ELSE  REWRITE(f); 

WRITELN(f,  JO); 

CLOSE(0; 

END; 

PROCEDURE  pdf_data_type.read_JO_from_file; 
VAR  f : TEXT; 

BEGIN 

ASSIGN(f,  ’JO.dat’); 

RESET(f); 

WHILE  NOT  EOF(f) 

DO  READLN(f,  JO); 

CLOSE(f); 

END; 

PROCEDURE  pdf_data_type.read_pdf_file; 
CONST  m_j)er_um  = le-6; 

VAR  f : TEXT; 
line  : STRING[80]; 
n,  check,  how_many  : INTEGER; 

X : DOUBLE; 
pdf_file_path  : PATHSTR; 
dir  : DIRSTR; 
name  ; NAMESTR; 
ext  : EXTSTR; 

BEGIN 

FSPLIT(model_prm_file_path,  dir,  name,  ext); 
pdf_file_path  :=  dir  + ’temp.pdf; 

ASSIGN(f,  pdf_file_path); 

RESET(f); 

FOR  n ;=  1 TO  8 
DO  READLN(f,  line); 

READLN(f,  how_many); 

FOR  n :=  1 TO  how_many 
DO  READLN(f,  line); 

READLN(f,  check); 

FOR  n :=  1 TO  how_many 
DO  READLN(f,  line); 

READLN(f,  check); 

X :=  -1; 
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WHILE  X < depth *m_per_um 
DO  IF  Vr  = 0 
THEN  BEGIN 
JO  :=  0; 

READLN(f,  X,  Jt); 

END 

ELSE  IF  shutter  = ’closed’ 
THEN  READLN(f,  x,  JO) 
ELSE  READLN(f,  x,  Jt); 
CLOSE(f); 

IF  Vr  <>  0 

THEN  IF  shutter  = ’closed’ 
THEN  write_J0_to_file 
ELSE  read_JO_from_file; 
END; 


PROCEDURE  pdf_data_type.append_QE_file; 

CONST  lambda_0  = 1239.85;  {Lambda  in  nm  at  which  1 W gives  1 A at  100%  QE.} 
VAR  f : TEXT; 

Jnet,  JlOO  : DOUBLE; 

BEGIN 

Jnet  :=  ABS(Jt-JO); 

JlOO  :=  lambda* phi/lambda_0; 

QE  :=  Jnet/JlOO; 

ASSIGN(f,  QE_file_path); 

APPEND(f); 

CASE  key  OF 

’O’  : WRITELN(f,  Vo:9:6,  ’ QE:9:6); 

’R’  : WRITELN(f,  Vr:9:6,  ’ QE:9:6); 

’W’  : WRITELN(f,  lambda:9:3,  ’ QE:9:6); 

’P’  : WRITELN(f,  energy:9:6,  ’ QE:9:6); 

’L’  : WRITELN(f,  phi:12,  ’ QE:9:6); 

END;  {CASE} 

CLOSE(f); 

END; 


BEGIN  {main 


WITH  pdf_data 
DO  BEGIN 
initialize; 

WRITELN(’READ_PDF.EXE  is  reading  TEMP.PDF.’); 
read_pdf_file; 

IF  shutter  = ’open’ 

THEN  BEGIN 

WRITELN(’READ_PDF.EXE  is  writing  ’,  QE_file_path); 
append_QE_file; 

END; 
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END; 

END. 
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APPENDIX  III.  A sample  of  the  screen  produced  by  MAXE_PRxM.EXE  when 
RUN_PC1D.BAT  is  executed. 

Parameter  file  (including  PRM  extension)  defining  photodiode: 

1337_new.prm 

Complete  path  to  file  in  which  to  store  QE  results: 
test.dat 

Type  of  experiment  to  simulate? 

(O) xide  bias  (Volts)  scan? 

(R)everse  bias  (Volts)  scan? 

(W)avelength  (nm)  scan? 

(P) hoton  energy  (eV)  scan? 

(L)inearity  (W/cm^2)  scan? 

Enter  letter  here:  w 

Start,  stop,  and  stop  values  from  (K)eyboard  or  all  values  from  (F)ile?  k 

First  value  for  scanned  parameter:  400 

Stop  value  for  scanned  parameter:  900 

Step  size  for  scanned  parameter:  10 

Incident  irradiance  (W  per  sq.  cm):  le-5 

Oxide  bias  voltage  (V):  0 

Reverse  bias  voltage  (V):  0 

Non-zero  thickness  of  front  surface  oxide  (nm):  25 

Depth  at  which  to  get  total  photocurrent  (um):  100 


33 


Part  II:  INTERPRETING  OXIDE-BIAS  EXPERIMENTS 


Jon  Geist 

NIST,  Gaithersburg,  MD  20899 
Rainer  Kohler  and  Roland  Goebel 
BIPM,  F-92310  Sevre,  France 
A.  M.  Robinson  and  C.  R.  James 
University  of  Alberta,  Alberta,  Canada  T6G  2G7 

Abstract 

The  semiconductor  device  modeling  program  PC-lD  and  the  programs  that  support  its 
use  in  high-accuracy  modeling  of  photodiodes,  ail  of  which  were  described  in  Part  I of 
this  paper,  are  used  to  simulate  oxide-bias  self-cedibration  experiments  on  three  different 
types  of  silicon  photodiodes.  It  is  shown  that  these  simulations  can  be  used  to  determine 
photodiode  characteristics,  including  the  internal  quantum  efficiency  for  the  different  types 
of  photodiodes.  In  the  latter  case,  the  simulations  provide  more  accurate  values  than  can 
be  determined  by  using  the  conventional  data  reduction  procedure,  and  an  uncertainty 
estimate  can  be  derived.  Finally,  it  is  shown  that  0.9997  ± 0.0003  is  a nominal  value  for 
the  internal  quantum  efficiency  of  one  type  of  photodiode  over  the  440  to  460  nm  spectral 
region. 

Key  words:  induced  junction  photodiode;  inversion  layer  photodiode;  numerical  modeling; 
numerical  simulation;  oxide-bias  experiment;  quajitum  efficiency 

1.  Introduction 

Part  II  presents  the  results  of  high- accuracy  simulations  of  various  oxide-bias  experiments 
performed  on  silicon  photodiodes.  The  simulations  were  carried  out  by  using  the  semicon- 
ductor device  modeling  program  PC-lD  and  a set  of  three  programs  designed  to  support 
the  use  of  PC-lD  in  this  specific  application.  Both  PC-lD  and  the  support  programs  are 
described  in  Part  I. 

The  oxide-bias  experiment  [1-3]  consists  of  recording  the  ratio  of  the  photocurrent  as 
a function  of  oxide-bias  voltage  to  the  zero-bias  photocurrent  during  irradiation  of  the 
photodiode  with  a stable  photon  fiux.  The  o^dde  bias  is  usually  applied  with  a transparent 
water-drop  electrode.  A review  is  given  in  Ref.  [3]. 

Part  II  is  organized  as  follows:  Section  2 describes  the  simulation  of  an  oxide-bias,  self- 
calibration experiment  on  an  EG&G  UV444B  type  photodiode.  The  purpose  of  this  sim- 
ulation is  to  demonstrate  the  derivation  of  an  accurate  value  for  the  zero-bias  quantum 
efficiency  including  an  uncertainty  estimate.  Section  3 describes  the  simulation  of  an  oxide- 
bias  experiment  on  a Hamamatsu  1337  type  photodiode.  The  purpose  of  this  simulation  is 
to  determine  nominal  values  for  the  charge  trapped  in  the  oxide  covering  the  front  surface 
of  the  photodiode,  and  for  the  surface  recombination  velocity  at  the  oxide-silicon  interface. 
Section  4 describes  the  simulation  of  an  oxide-bias  experiment  on  a UDT  UVlOO  type 
photodiode.  The  purpose  of  this  simulation  is  to  determine  a nominal  high-accuracy  value 
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for  the  internal  quantum  efficiency  of  this  type  of  photodiode  and  the  nominal  wavelength 
region  where  this  value  is  valid.  Section  5 gives  the  conclusions  of  Part  II. 

2.  EG&:G  UV444B  Photodiode 

The  EG(kG  UV444B  silicon  photodiode  was  the  first  of  its  type  to  be  self-calibrated  by 
using  the  oxide-bias  experiment  [2] . This  type  of  photodiode  consists  of  a 300  ^m  thick,  300 
Q.  cm,  phosphorus-doped,  n-type  silicon  substrate  with  an  n^-type  phosphorus  diffusion 
in  the  rear  surface,  a p^-type  boron  diffusion  in  the  front  surface,  and  a nominally  120  nm 
thick,  thermedly  grown  oxide-passivation  layer  over  the  front  surface.  Figure  1 shows  the 
front-region  doping  profile  that  was  obtained  from  spreading  resistance  measurements  on  a 
representative  UV444B  photodiode.  Figure  1 also  shows  the  equilibrium  hole  concentration 
calculated  for  this  doping  distribution  by  using  PC-lD. 

The  rear-region,  n^-type  diffusion  in  the  UV444B  was  not  modeled;  instead  it  was  assumed 
that  the  recombination  velocity  at  the  rear  of  the  device  was  zero.  This  assumption  is  not 
physically  realizable,  but  it  results  in  a mirror  for  minority  carriers  at  the  rear  of  the 
device,  just  as  does  the  n'^-type  diffusion.  This  assumption  is  useful  because  it  requires 
fewer  finite  elements,  thereby  leaving  more  of  them  available  for  modeling  the  front  region 
and  depletion  region. 

Figure  2 compares  simulated  and  experimental  results  of  an  oxide-bias  experiment  on  a 
UV444B  photodiode  for  irradiation  by  a 3 mm  diameter  spot  of  476.2-nm  laser  radiation 
through  a water-drop  electrode.  Reference  [4]  describes  the  apparatus  used  to  obtain  the 
experimental  data.  The  absorption-coefficient  data  in  the  file  SIL.WEAK.ABS  used  for 
this  simulation  were  calculated  from  an  equation  [5]  that  was  fit  to  the  data  measured  by 
Weakliem  and  Redfield  [6].  These  data  are  described  in  Part  I.  Both  the  surface  recombi- 
nation velocity  S and  the  number  density  N„  of  the  oxide  fixed  charge  were  adjusted  in 
the  simulation  to  obtain  a good  fit  to  the  experimental  data  in  figure  2. 

The  procedure  used  to  fit  the  simulated  data  to  the  experimental  data  shown  in  figure  2 
was  as  follows:  first,  a value  of  S was  chosen,  then  the  value  of  N,,  was  adjusted  until 
the  average  residual  of  the  fit  to  the  photocurrent  ratio  was  zero  for  the  last  two  data 
points  in  the  figure.  This  gave  a value  of  N„  for  each  value  of  S that  was  chosen.  The 
best  combination  of  S and  N,,  was  then  obteiined  by  minimizing  the  residual  standard 
deviation  over  all  of  the  data  shown  in  figure  2. 

Figures  3,  4,  and  5 plot  the  residuals  (differences)  between  the  results  of  the  fits  and  the 
experimental  results  and  table  1 lists  the  residual  standard  deviation  for  three  different 
choices  of  S and  N,j.  It  is  clear  by  inspection  of  the  figures  that  the  combination  of  S and 
N,,  used  for  figure  3 gives  a better  overall  fit  than  either  of  the  other  choices,  even  though 
the  residual  standard  deviation  for  figure  3 is  not  very  much  lower  than  that  for  figure  4 or 
figure  5.  One  reason  for  the  relatively  small  decrease  in  residual  standard  deviation  with 
the  best  fit  is  the  effect  of  the  apparent  outlier  points  near  -1  V and  -10.5  V oxide  bias. 

The  assignment  of  S = 6000  cm/s  and  N„  = 1.942x10^^  cm“^  with  estimated  uncertainties 
of  ±1000  cm/s  and  ±0.070  x 10^^  cm“^,  respectively,  is  justified  by  figures  3,  4,  and  5.  This 
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Figure  1.  The  front-region  boron  concentration  (continuous  line)  and  the  equilibrium 
hole  concentration  (dashed  line)  calculated  by  PC-lD  for  a typical  EG<SiG  UV444B  type 
photodiode. 
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Figure  2.  Ratio  of  the  photocurrent  I(Vo)  as  a function  of  oxide-bias  voltage  Vq  to  the 
photocurrent  1(0)  at  zero-bias  voltage  as  measured  on  an  EG&G  UV444B  silicon  photo- 
diode and  as  calculated  with  PC-lD  with  a front-surface  recombination  velocity  S = 6000 
cm/s  and  an  oxide  trapped  charge  number  density  Nj,  = 1.942  x 10^^  cm“^. 
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RESIDUAL 


Figure  3.  The  residuals  between  the  experimental  data  of  figure  2 and  the  simulated  data 
calculated  by  PC-lD  when  S = 6000  cm/s  and  N,,  = 1.942  x 10^^  cm~^. 
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RESIDUAL 


Figure  4.  The  residueds  between  the  experimental  data  of  figure  2 and  the  simulated  data 
calculated  by  PC-lD  when  S = 5000  cm/s  and  N,,  = 2.011  x 10^^  cm“^. 
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RESIDUAL 


Figure  5.  The  residuals  between  the  experimental  data  of  figure  2 and  the  simulated  data 
calculated  by  PC-lD  when  S = 7000  cm/s  and  N,,  = 1.883  x 10^^  cm“^. 
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Table  1.  Comparison  of  simulated  oxide  bias  results  with  experimental  results  for  a 
UV444B  type  photodiode  where  Res.  Std.  Dev.  is  the  standard  deviation  of  the  dif- 
ference between  simulated  and  experimental  results. 


S (cm/s) 

N,,  (10^2  cm 

Res.  Std.  D 

5000 

2.011 

0.002233 

6000 

1.942 

0.002076 

7000 

1.883 

0.002093 

41 


assignment  corresponds  to  an  internal  quantum  deficiency  (one  minus  the  internal  quantum 
efRciency)  at  zero  oxide  bias  of  0.1690  ±0.0021.  The  stated  uncertainty  is  the  maximum 
difference  between  the  internal  quantum  deficiency  calculated  for  the  combination  of  S 
and  N,,  used  in  figure  3,  and  those  calculated  for  the  combinations  of  S and  N,,  used  for 
figures  4 and  5. 

The  assumption  that  the  average  of  the  last  two  data  points  in  figure  3 corresponds  to  1009c 
internal  quantum  efficiency  gives  a zero-bias  internal  quantum  deficiency  of  0.1634  with 
no  reliable  way  to  estimate  the  uncertainty  associated  with  that  value.  This  assumption, 
which  is  usually  used  with  the  oxide-bias,  self- calibration  experiment,  produces  a result  that 
differs  from  the  more  accurate  value  determined  from  the  fit  by  0.0056.  This  difference 
is  about  3.3%  of  the  effect  being  calculated,  and  more  than  2.5  times  the  uncertainty 
associated  with  the  more  accurate  value. 

It  must  be  emphasized  that  the  results  reported  here  are  illustrative  rather  than  definitive 
since  the  doping  distribution  data  in  figure  1 were  obtained  from  one  photodiode,  and 
the  oxide-bias  data  in  figure  2 were  obtained  from  a different  photodiode.  An  additional 
uncertainty  associated  with  the  effect  of  variations  in  doping  distribution  would  be  needed, 
or  the  doping  distribution  must  be  measured  on  the  actual  photodiode  used  in  the  self- 
calibration. 

Because  oxide  bias  voltages  large  enough  to  nearly  saturate  the  photodiode  response  are 
also  large  enough  to  change  S or  N,,,  the  oxide-bias  measurement  is  already  destructive 
to  a certain  extent,  and  it  must  be  carried  out  as  described  in  Ref.  [7]  in  order  to  yield 
the  highest  accuracy  of  which  it  is  capable.  Therefore,  the  sacrifice  of  the  photodiode 
for  a spreading  resistance  measurement  is  not  as  impractical  as  it  might  at  first  seem. 
Nevertheless,  it  is  inconvenient,  and  more  convenient  approaches  are  being  developed  [3;. 

The  results  reported  here  might  be  considered  illustrative  for  another  reason.  The  fitting 
and  uncertainty  analysis  were  carried  out  by  hand  in  a somewhat  subjective  manner. 
Ideally,  the  fitting  and  uncertainty  analysis  would  be  carried  out  automatically  by  a fitting 
routine  incorporated  into  the  set  of  programs  used  to  run  PC-lD  for  photodiode  modeling. 
However,  this  task  might  not  be  as  straightforward  as  it  first  would  seem  because  the  fit 
should  be  constrained  to  pass  through  the  mean  of  the  last  few  data  points  at  the  highest 
oxide  bias  voltage,  as  it  was  in  the  hand  fit  reported  here.  Without  this  constraint,  the 
minimum  residual  standard  deviation  might  be  obtained  for  a choice  of  parameters  that 
clearly  does  not  fit  the  data  well  at  the  highest  oxide  bias  voltages.  It  is  clear  that  the 
data  at  the  highest  oxide  bias  voltages  should  be  weighted  more  heavily  in  the  fit,  but  it 
is  not  clear  how  much  more  heavily.  Therefore,  an  automated  program  does  not  eliminate 
the  subjectivity,  but  only  transfers  it  to  the  weighting  of  the  data. 

3.  Hamamatsu  1337  Photodiode 

Like  the  EG&G  UV444B  photodiode,  the  Hamamatsu  1337  is  a p^nn"^  photodiode,  but 
its  front-region  doping  profile  is  shallower,  its  background  n-type  dopant  concentration  is 
5 X 10^^  cm“^,  and  its  oxide-passivation  layer  has  a nominal  thickness  of  25  nm.  Figure  6 
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Figure  6.  The  relative  front -region  boron  concentration  for  six  typical  Hamamatsu  1337 
type  photodiodes.  The  two  profiles  shown  with  solid  lines  bracket  the  other  four  in  their 
effect  on  the  shape  of  the  internal  quantum-deficiency  spectrum. 
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shows  relative  front-region  doping  profiles  for  six  different  Hamamatsu  1337  photodiodes 
taken  from  four  different  batches.  The  curves  in  figure  6 were  normalized  to  unity  at 
the  oxide-silicon  interface  because  it  is  the  shape  of  the  doping  profile  rather  than  the 
magnitude  that  has  the  major  effect  on  the  internal  quantum  efficiency  [8].  The  measured 
dopant  concentrations  at  the  oxide-silicon  interface  varied  between  2 x 10^®  and  1 x 10^^ 
cm~^  for  the  six  photodiodes  of  figure  6. 

Figure  7 shows  the  results  of  modeling  an  oxide-bias  experiment  on  a 1337  type  photodiode 
at  476.2  nm  by  using  the  doping  profile  that  corresponds  to  the  lower  solid  curve  in 
figure  6.  The  fit  shown  in  figure  7 was  obtained  with  N,,  = —3.0  x 10^^  cm~^  and  S = 
1.835  X 10^  cm/s,  and  gave  a zero-bieis  internal  quantum  deficiency  of  0.003199.  A fit  that 
looked  virtually  identical  to  that  shown  in  figure  7 was  obtained  for  the  doping  profile 
corresponding  to  the  upper  solid  curve  in  figure  6 by  using  N,,  = —4  x 10^^  cm“^  and  S 
= 1.56  X 10^  cm/s.  This  fit  gave  a zero- bias  internal  quantum  efficiency  of  0.003280.  As 
in  the  case  of  the  UV444B  type  photodiode,  the  SIL.WEAK.ABS  absorption-coefficient 
data  were  used,  and  the  rear  of  the  photodiode  was  modeled  as  uniformly  doped  with  zero 
recombination  velocity. 

Figure  8 compares  the  internal  quantum  deficiencies  calculated  from  the  doping  profile 
corresponding  to  the  solid  curves  in  figure  6 by  using  the  values  of  N,,  and  S that  were 
obtained  from  the  simulated  oxide-bias  experiments  based  on  each  of  those  profiles.  The 
fact  that  the  two  curves  have  very  similar  shapes  is  very  important.  It  suggests  that  a 
single  function  of  wavelength  might  be  used  to  extrapolate  the  internal  quantum  efficiency 
of  any  1337  type  photodiode,  rather  than  a different  function  for  each  different  photodiode. 
This  possibility  is  examined  in  detail  in  Part  III. 

Notice  that  the  sign  of  the  trapped  charge  in  the  1337  type  photodiode  is  the  opposite  of 
that  for  the  UV444B  photodiode,  and  that  the  interface  recombination  velocity  is  much 
larger.  Both  the  sign  of  the  trapped  charge  and  the  large  value  of  the  recombination 
velocity  are  consistent  with  an  oxide  prepared  by  chemical  vapor  deposition  rather  than 
by  thermal  oxidation,  because  the  former  is  known  to  produce  trapped  positive  charge  and 
a low  surface  recombination  velocity.  The  modeling  carried  out  here  shows  that  a large 
negative  charge  in  the  oxide  more  than  makes  up  for  the  greater  surface  recombination 
velocity  by  repelling  minority  carriers  from  the  surface. 

The  fact  that  the  1337  type  photodiode  has  negative  rather  than  positive  charge  trapped 
in  its  oxide  is  important  for  its  possible  use  in  quantum-efficiency  interpolation  and  ex- 
trapolation. If  the  oxide  contained  a significant  amount  of  positive  charge,  as  assumed 
in  Ref.  [9],  the  oxide-silicon  interface  would  be  depleted  of  holes,  leading  to  a very  large 
decrease  in  the  energy  needed  to  impact-ionize  a valence  band  electron.  However,  since 
the  oxide  contains  negative  charge,  the  interface  is  accumulated.  In  this  case,  even  though 
the  charge  density  is  quite  large,  there  is  little  or  no  decrease  in  the  energy  needed  to 
impact-ionize  a valence  band  electron  [9].  This  is  important  because  a large  decrease  in 
the  energy  required  to  impact-ionize  a valence  band  electron  might  cause  the  quantum 
yield  [10]  for  electron-hole  pair  production  of  the  silicon  near  the  front  surface  to  exceed 
unity  at  wavelengths  as  long  as  450  nm.  This  would  limit  the  usefulness  of  this  type 
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Figure  7.  Ratio  of  the  photocurrent  I(Vo)  as  a function  of  oxide-bias  voltage  Vq  to  the 
photocurrent  1(0)  at  zero-bias  voltage  as  measured  on  a Hamamatsu  1337  silicon  photodi- 
ode and  as  calculated  with  PC-lD  for  the  doping  distribution  corresponding  to  the  upper 
solid  curve  in  figure  6 with  S = 1.835  cm/s  and  N,,  = —3.0  x 10^^  cm“^. 
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Figure  8.  Simulated  internal  quantum  deficiencies  for  the  1337  type  photodiode  having 
the  oxide-bias  data  shown  in  figure  7:  (solid  curve)  for  the  doping  profile  that  corresponds 
to  the  lower  solid  curve  in  figure  6 with  Nj,  = —3.0  x 10^^  and  S = 1.835  x 10®  cm/s,  and 
(dashed  curve)  for  the  doping  profile  that  corresponds  to  the  upper  solid  curve  in  figure  6 
with  N,,  — —4  X 10^^  and  S = 1.56  x 10®  cm/s. 
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of  photodiode  for  extrapolating  and  interpolating  quantum-efficiency  calibrations.  If  this 
limitation  were  encountered,  it  could  only  be  overcome  with  a very  accurate  model  of  the 
quantum  yield  in  the  presence  of  surface  fields.  Such  a model  appears  well  beyond  the 
current  state  of  the  art  [9]. 

Even  though  the  uncertainties  in  the  internal  quantum  deficiency  associated  with  the 
diode-to-diode  variations  in  doping  profile  are  quite  small  for  the  1337  type  photodiode, 
the  oxide-bias  experiment  is  still  less  than  ideal  for  self-calibration  [3].  However,  since  the 
1337  is  so  well  suited  for  interpolating  and  extrapolating  internal  quantum-deficiency  data, 
as  shown  in  Part  III,  all  that  is  needed  are  ways  to  calibrate  the  photodiode  that  work 
over  limited  spectral  regions  near  450  nm  and  near  900  nm.  The  next  section  illustrates 
how  simulations  can  be  used  to  predict  the  internal  quantum  efficiency  of  a different  type 
of  photodiode  with  an  uncertainty  of  about  0.0003  at  450  nm. 

4.  UDT  UVlOO  Photodiode 

The  UDT  UVlOO  type  photodiode  is  usually  assumed  to  have  an  internal  quantum  effi- 
ciency of  unity  in  the  short  wavelength  portion  of  the  visible  [11].  In  fact,  this  assumption 
is  the  basis  for  the  use  of  the  QED  100  and  QED  200  multiple  reflection  radiometers 
[12],  which  contain  modified  UVlOO  photodiodes,  as  high-accuracy  absolute  radiometric 
standards. 

As  the  state  of  the  art  in  absolute  radiometry  improves  with  time,  the  assumption  of  unity 
quantum  efficiency  for  the  UVlOO  type  photodiode  becomes  more  questionable.  What 
is  needed  is  a nomined  internal  quantum-deficiency  spectrum  accompanied  by  a reliable 
uncertainty  spectrum.  Since  the  1337  type  photodiode  can  be  used  to  interpolate  or  ex- 
trapolate internal  quantum-efificiency  calibrations,  as  shown  in  Part  III,  it  is  only  necessary 
to  have  these  spectra  over  a limited  spectral  region  above  400  nm.  (Below  400  nm,  the 
non-unity  quantum  yield  would  be  a problem.)  One  way  to  obtain  these  spectra  is  to  cal- 
culate them  from  nominal  characteristics  of  the  UVlOO  type  of  photodiode.  The  pertinent 
characteristics  can  be  obtained  by  fitting  numerical  simulations  of  oxide-bias  measurements 
to  experimental  results  as  demonstrated  in  the  preceding  sections. 

The  UVlOO  type  of  photodiode  is  quite  different  from  the  UV444B  and  1337  type  photodi- 
odes. Its  n^p  junction  is  not  formed  by  introducing  a front-region  dopant,  but  is  induced 
in  a 100  H cm  p-type  (approximately  1.3  x 10^**  cm“^,  boron-doped)  silicon  substrate  by 
the  growth  of  an  oxide  passivation  layer  having  a large  concentration  of  trapped  positive 
charge  [13]. 

Figure  9 compares  the  results  of  simulated  and  experimental  oxide-bias  experiments  on  a 
UVlOO  type  photodiode.  Here  again  the  absorption-coefficient  data  in  SIL.WEAK.ABS 
were  used,  and  the  rear  of  the  photodiode  was  modeled  as  uniformly  doped  with  a surface 
recombination  velocity  of  zero.  Since  the  UVlOO  type  photodiode  is  an  n'^p  type  rather 
than  a p'^’n  type,  the  effect  of  increasing  negative  bias  is  to  decrease  rather  than  increase 
the  photocurrent.  However,  because  the  zero-bias  quantum  efficiency  is  so  close  to  unity, 
this  procedure  samples  the  portion  of  the  curve  with  the  most  information  about  the  values 
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I (Vo)/ 1(0) 


Figure  9.  Ratio  of  the  photocurrent  I(Vo)  as  a function  of  oxide-bias  voltage  Vq  to  the 
photocurrent  1(0)  at  zero- bias  voltage  as  measured  on  a typicad  UDT  UVlOO  silicon  pho- 
todiode and  as  calculated  with  PC-lD  for  S = 3.0  x 10*  cm/s  and  N,,  = 1.38  x 10^^ 
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of  N,,  and  S,  as  can  be  seen  in  the  figure. 

These  values  of  N,,  and  S are  used  to  simulate  the  internal  quantum-deficiency  spectra  for 
a nominal  UVlOO  photodiode  in  figure  10.  The  three  different  curves  shown  in  figure  10 
correspond  to  minority  carrier  lifetimes  r of  1 ^s,  10  /xs,  and  81  ms,  respectively,  in  the  rear 
region  of  the  photodiode.  The  lifetimes  of  1 fis  and  10  /xs  give  internal  quantum  deficiencies 
at  700  nm  that  bracket  the  range  measured  for  six  typical  UVlOO  type  photodiodes  T4\ 
and  the  81  ms  lifetime  is  the  longest  lifetime  that  the  default  silicon  material  model  of 
PC-lD  will  allow  for  100  Cl  cm,  p-type  silicon  with  a recombination  trap  at  mid-gap. 

The  region  between  the  curves  corresponding  to  r = 1 and  r = 10  /xs  in  figure  10  is  a 
nominal  internal  quantum-deficiency  spectrum  that  applies  to  UVlOO  type  photodiodes 
and  to  QED  100  and  QED  200  radiometers  at  zero  reverse  bias  [3,12].  Figure  11  compares 
the  spectrum  for  r = 1 /xs  around  its  minimum  in  the  400  to  500  nm  spectral  region 
with  spectra  obtained  for  the  cases  where  either  the  front-region  surface  recombination 
velocity  is  doubled,  or  the  positive  charge  trapped  in  the  oxide  is  halved.  These  are 
conservative  limits  to  associate  with  the  mass  production  of  photodiodes  in  a controlled 
planar-silicon  fabrication  process.  These  should  also  be  conservative  limits  for  the  effects 
of  environmental  stress,  such  as  humidity,  temperature,  and  irradiation  with  ultraviolet 
radiation,  for  UVlOO  photodiodes  [15-16]  made  after  1986. 

The  curves  in  figure  11  can  be  used  to  estimate  uncertainties  to  be  associated  with  the 
nominal  values  of  the  internal  quantum  efficiency  of  the  UVlOO  type  photodiode.  Summing 
in  quadrature  the  differences  between  the  nominal  spectrum  and  the  other  two  spectra  in 
figure  11  produces  a somewhat  conservative  estimate  since  increases  in  surface  recombina- 
tion velocity  tend  to  be  correlated  with  increases  in  trapped  charge.  In  the  region  from 
440  to  460  nm,  a nominal  value  of  0.0003  ± 0.0003  is  a practical  summary  of  the  results 
of  such  an  analysis.  This  result  should  be  considered  a one-standard-deviation  equivalent, 
rather  than  a limit  of  error,  since  the  whole  analysis  is  based  on  the  results  obtained  for  a 
single  device.  Outside  the  440  to  460  nm  spectral  region  the  uncertainties  grow  so  rapidly 
that  they  axe  not  really  practical  for  high-accuracy  applications. 

There  are  a number  of  sources  of  error  not  yet  considered,  but  since  their  effect  scales  with 
the  magnitude  of  the  internal  quantum  deficiency,  they  are  edl  negligible  with  respect  to 
±0.0003.  This  is  illustrated  in  figure  12  for  the  uncertainty  associated  with  the  absorption 
coefficient  of  silicon  in  the  440-to  460-nm  spectral  region,  which  is  the  largest  uncertainty 
not  considered  so  far. 

Figure  12  compares  the  spectrum  of  figure  11  for  r = 1 /xs,  which  was  based  on  the 
absorption-coefficient  data  in  SIL_WEAK.ABS,  with  the  same  spectra  calculated  for  iden- 
tical parameters,  but  based  on  absorption-coefficient  data  (measured  by  Philipp  [17],  and 
described  in  Paxt  I)  that  is  stored  in  SIL_PHIL.ABS.  When  added  in  quadrature  to  the 
±0.0003  assigned  to  the  variations  in  N„  and  S from  diode  to  diode,  the  effect  of  the  un- 
certainty in  the  absorption  coefficient  is  negligible.  This  might  not  remain  true  if  a much 
larger  wavelength  interval  were  considered. 
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Figure  10.  Internal  quantum-deficiency  spectra  calculated  by  PC-lD  for  S = 3.0  x 10^ 
cm/s,  N,,  = 1.38  X 10^^  cm“^,  and  rear-region  lifetimes  r of  81  ms  (solid  line),  10  /is 
(dashed  line),  and  1 /xs  (dot-dashed  line). 
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Figure  11.  Internal  quantum-deficiency  spectra  calculated  by  PC-lD  for  a rear-region 
lifetime  r = 1 /is  with  S = 3.0  x 10®  cm/s  and  N,,  = 1.38  x 10^^  cm”^  (solid  line),  with 
r = 1 |xs  with  S = 6.0  x 10®  cm/s  and  N,,  = 1.38  x 10^^  cm“^  (dashed  line),  and  whth  r 
= 1 /zs  with  S = 3.0  X 10®  cm/s  and  N,,  = 0.69  x 10^^  cm“^  (dot-dashed  line). 
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Figure  12.  Comparison  of  the  internal  quantum-deficiency  spectra  calculated  for  a UVlOO 
type  photodiode  for  a rear-region  lifetime  r = 1 ^s  with  S = 3.0  x 10^  cm/s  and  N,,  = 
1.38  X 10^^  cm~^  with  the  absorption- coefficient  data  in  the  SIL.WEAK.ABS  file  (solid 
line),  and  with  the  absorption-coefficient  data  in  the  SIL_PHIL.ABS  file  (dashed  line).  The 
data  in  these  files  are  described  in  more  detail  in  Part  I. 
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Because  of  the  large  long- wavelength  quantum  deficiency  shown  in  figure  10,  reverse  bias 
is  often  applied  to  UVlOO  photodiodes  and  QED  100  and  200  radiometers  vrhen  they  are 
used  at  the  longer  wavelengths.  By  moving  the  depletion  region  toward  the  rear  of  the 
photodiode,  the  reverse  bias  forces  the  internal  quantum  deficiency  to  approach  the  lower 
quantum  deficiency  limit  given  by  the  81  ms  curve  in  figure  10  [12,18-19].  Reverse  bias 
is  also  applied  to  eliminate  the  saturation  type  nonlinearity  that  is  sometimes  observed 
with  photocurrents  above  a few  microamperes.  However,  recent  results  14]  suggest  that 
moderate  bias  levels  can  cause  gain,  resulting  in  an  internal  quantum  efficiency  in  excess  of 
unity.  This  effectively  precludes  the  application  of  moderate  reverse  bias  in  high-accuracy 
applications. 

Figure  13  plots  a simulation  of  the  decrease  in  the  internal  quantum  deficiency  as  a function 
of  reverse  bias  for  the  UVlOO  type  photodiode  at  the  442,  476,  and  514  nm  laser  lines. 
The  feature  in  the  curves  for  reverse-bias  voltages  less  (more  negative)  than  -3.5  V is  a 
symptom  of  a problem  with  PC-lD.  It  may  mean  that  more  than  150  finite  elements  are 
needed  to  maintain  the  level  of  accuracy  achieved  as  the  reverse  bias  is  increased,  but  this 
is  not  certain.  In  any  case,  it  not  only  invalidates  the  data  for  reverse-bias  voltages  more 
negative  than  -3.5  V,  but  it  also  casts  some  doubt  on  any  conclusions  drawn  about  the 
data  at  the  less  negative  bias  voltages  as  well.  At  best,  figure  13  suggests  that  480  nm  is  a 
practical  upper  limit  to  the  wavelength  at  which  reverse  bias  can  be  used  to  eliminate  the 
effect  of  recombination  in  the  rear  region  of  the  photodiode  without  risking  gsun.  However, 
the  best  advice  would  be  never  to  use  a UVlOO  photodiode  or  QED  radiometer  with  a bias 
voltage  without  first  assuring  experimentally  that  gain  is  not  a problem. 

5.  Conclusion 

Part  II  of  this  paper  has  demonstrated  the  use  of  Version  2 of  PC-lD.EXE,  RUN_PClD.BAT 
MAKE_PRM.EXE,  and  READ_PDF.EXE  to  simulate  various  oxide-bias  experiments  of 
interest  to  high-accuracy  applications  of  silicon  photodiodes.  In  each  case,  the  simulation 
was  used  to  interpret  an  oxide-bias  experiment,  with  different  sorts  of  information  being 
sought  from  the  different  experiments. 

It  was  shown  that  the  simulations  can  be  used  to  derive  a more  accurate  internal  quantum- 
deficiency  value  from  an  oxide-bias  experiment  than  that  available  from  the  conventional 
data  reduction,  and  that  an  uncertainty  can  be  associated  with  the  value  derived  from  the 
simulation  based  on  how  well  the  simulated  oxide-bias  data  can  be  fitted  to  the  experi- 
mental data.  It  was  also  shown  that  simulations  can  be  used  to  determine  nominal  values 
for  the  front-surface  recombination  velocity  and  for  the  charge  trapped  in  the  front-surface 
oxide,  and  that  values  for  these  parameters  can  be  used  to  determine  nominal  internal 
quantum-deficiency  spectra  for  various  types  of  silicon  photodiodes. 

As  well  as  illustrating  these  general  ideas,  some  more  specific  results  were  also  obtained. 
Firstly,  it  was  shown  that  a nominal  internal  quantum  deficiency  of  0.0003  ± 0.0003  is 
appropriate  for  the  unbiased  UDT  UVlOO  type  photodiode  in  the  440  to  460  nm  spec- 
tral region.  Secondly,  it  was  shown  that  the  shape  of  the  internal  quantum-deficiency 
spectrum  of  the  1337  type  photodiode  is  quite  independent  of  typical  variations  in  doping 
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Figure  13.  Simulation  of  the  decrease  in  internal  quantum- deficiency  of  a UVlOO  type 
photodiode  with  reverse  bieis  at  442  nm  (heavy  solid  line),  468  nm  (light  solid  line),  and 
514  nm  (dashed  line).  The  values  of  N„,  S,  and  r are  as  in  figure  12.  The  feature  near 
-4  V is  an  artifact  of  the  simulation.  The  open  circle  and  closed  circle  are  the  internal 
quantum  deficiencies  at  442  nm  and  468  nm,  respectively,  for  r = 81  ms.  The  internal 
quantum  deficiency  at  514  nm  for  the  same  condition  is  too  small  to  appear  on  the  graph. 
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profile.  This  result  suggests  that  the  1337  type  photodiode  might  serve  to  extrapolate 
high-accuracy  calibrations  obtained  from  cryogenic  absolute  radiometers  or  from  QED 
type  radiometers  in  the  440  to  460  nm  spectral  region  to  longer  wavelengths  with  little  or 
no  loss  of  accuracy.  Part  III  of  this  paper  shows  that  this  is  the  case.  Finally,  it  was  shown 
that  reverse-bias  simulations  based  on  the  programs  described  in  Part  I of  this  paper  do 
not  work  well  enough  with  inversion  layer  (induced  junction)  photodiodes  to  be  useful  for 
high-accuracy  applications. 
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Abstract 

The  semiconductor  device  modeling  program  PC-lD  and  the  programs  that  support  its 
use  in  high-accuracy  modeling  of  photodiodes,  all  of  which  were  described  in  Part  I,  are 
used  to  simulate  the  interpolation  of  high- accuracy  internal  quantum-efficiency  calibrations 
in  the  spectral  region  between  450  nm  and  850  nm.  Convenient  interpolation  formulae 
that  depend  only  upon  wavelength  are  derived.  Uncertainty  spectra  for  a number  of 
sources  of  error  are  also  derived.  The  formulae  are  normalized  to  experimental  internal- 
quantum  efficiency  calibrations  in  the  440  to  470  nm  spectral  region  and  at  860  nm  and 
are  used  to  interpolate  the  calibration  values  between  these  wavelengths.  The  results  of 
the  interpolations  are  compared  with  experimental  calibration  data  that  are  available  at 
a few  wavelengths  between  440  and  860  nm.  The  disagreement  between  the  interpolated 
and  measured  internal  quantum-efficiency  data  is  never  worse  than  0.0003. 

Key  words:  extrapolation;  internal  quantum-efficiency;  interpolation;  photodiode;  silicon; 
visible 

1.  Introduction 

Part  II  of  this  paper  suggested  that  the  Hamamatsu  1337  type  photodiode  might  be  con- 
veniently used  for  interpolating  or  extrapolating  high-accuracy  quantum-efficiencies  in  the 
400  to  900  nm  spectral  region.  The  idea  of  quantum-efficiency  extrapolation  was  implicit 
in  the  first  detailed  description  of  how  a silicon  photodiode  could  be  used  (in  this  case 
in  conjunction  with  a thermal  detector  having  a spectredly  fiat  responsivity)  to  realize  a 
scale  of  absolute  spectral  responsivity  [1],  and  interpolation  was  explicitly  considered  in 
Ref.  [2].  Beyond  these  two  publications,  however,  little  progress  has  been  reported  in  this 
area  until  recently. 

Hoyt  et  al.  [3]  (using  early  results  of  the  work  described  here)  used  three  Hamamatsu 
1337  type  photodiodes  to  extrapolate  a quantum-efficiency  calibration  obtained  at  441.6 
nm  with  a UDT  QED  100  radiometer  to  633.0  nm,  and  reported  an  average  difference 
of  0.04%  compared  to  independent  measurements  based  on  their  electrically  calibrated 
cryogenic  cavity  radiometer.  Even  though  0.04%  is  quite  good  by  conventional  radiometric 
standards,  it  was  a factor  of  two  larger  than  the  combined  uncertainty  estimate  for  the  two 
measurement  techniques  being  compared,  and  further  study  was  proposed  by  Hoyt  et  al. 
More  recently,  Zalewski  and  Hoyt  [4]  have  reported  a more  direct  test  of  the  accuracy  of 
the  use  of  the  Hamamatsu  1337  type  photodiode  to  extrapolate  quantum  efficiency  values. 
The  results  of  this  test  agreed  to  within  the  combined  estimated  uncertainty,  and  two 
errors  were  discovered  in  the  earlier  work.  When  these  errors  were  corrected,  agreement 
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within  the  estimated  uncertainty  of  the  earlier  intercomparison  was  obtained. 

The  purpose  of  Part  III  is  to  develop  formulae  and  uncertainty  estimates  for  interpo- 
lating and  extrapolating  internal  quantum-efficiency  calibrations  with  Hamamatsu  1337 
photodiodes.  It  will  be  shown  that  this  type  of  photodiode  has  characteristics  that  suit 
it  particularly  well  for  this  task  over  the  440  to  860  nm  spectral  region.  The  remainder 
of  Part  III  is  organized  as  follows:  Section  2 investigates  the  errors  associated  with  the 
use  of  relations  that  are  only  approximate  for  the  1337  type  photodiode,  but  which  sim- 
plify interpolation  and  extrapolation  of  the  internal  quantum  deficiency  (one  minus  the 
internal  quantum  efficiency)  for  this  type  of  photodiode.  Section  3 investigates  the  various 
sources  of  error  associated  with  the  actual  calculation  of  internal  quantum-deficiency  spec- 
tra using  the  semiconductor-device  modeling  program  PC-lD  and  the  support  programs 
described  in  Part  1.  Finally,  section  4 compares  the  results  obtained  here  with  experimental 
measurements  and  demonstrates  agreement  within  0.0003. 

2.  Errors  Associated  with  the  Quantum-Deficiency  Approximations 

For  the  purposes  of  extrapolating  and  interpolating  the  internal  quantum  efficiency  of 
Hamamatsu  1337  type  photodiodes,  the  internal  quantum  deficiency  can  be  approximated 
by  a function  S{X,  iV,,,  5,  r^)  of  wavelength  A and  of  three  recombination- related  variables; 
Ngt  the  charge  number  density  stored  in  the  front-surface  oxide,  5 the  surface  recombina- 
tion velocity  at  the  oxide-sihcon  interface,  and  the  minority-carrier  lifetime  in  the  rear 
region. 

For  silicon  photodiodes,  ^(A,  Ngg^  5,  Tr)  can  be  approximated  by  a function  that  has  proper- 
ties that  make  it  particularly  useful  for  extrapolating  and  interpolating  quantum  deficiency 
measurements.  The  approximate  equation  for  ^(A,  iV,,,  5,  r^)  is 


S{X,Ngg,S,Tr)  = Kf{N,g,S)Sf{X) 


+ Kr{Tr)6r{X), 


(1) 


where 


«r(A)  = 0, 


(2) 


for  A less  than  some  wavelength  Aq. 

To  see  how  this  is  useful  for  interpolating  quantum  deficiency  values,  suppose  that  6{X  f , Ng,, 
5,  Tr)  has  been  measured  at  wavelength  Xf  < Aq,  and  ^(Ar,  iV,,,  S,  r^)  has  been  measured 
at  wavelength  Ar  > Aq.  The  values  for  K f{Ntg,  S)  and  Kr{Tr)  can  then  be  determined 
from 
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(3a) 


K,{N..,S)  = 


h(h) 


and 


^ _ S{\r,N..,S,r,)-Kf{N,.,S)Sf{K) 

Sr{K) 

and  these  values  can  be  used  in  eq  (1)  to  interpolate  the  quantum  deficiency  between  Xf 
and  Xr.  In  the  above  equations,  Sf{X)  and  <^r(A)  are  functions  only  of  wavelength,  Xf  and 
A,,  are  the  wavelengths  at  which  the  internal  quantum  deficiency  is  to  be  normalized,  and 
Ao  will  generally  fall  between  600  nm  and  650  nm.  The  subscript  / stands  for  the  front 
of  the  photodiode  since  this  is  the  location  of  the  recombination  mechanisms  that  cause 
Sf{X)  to  exceed  zero.  The  subscript  r stands  for  the  rear  of  the  photodiode  since  this  is 
the  location  of  the  recombination  mechanism  that  cause  ^r(A)  to  exceed  zero.  Therefore, 
Sf{X)  and  ^r(A)  will  be  called  the  front-region  and  rear-region  internal  quantum-deficiency 
spectra,  respectively.  The  functional  dependence  oi  Kf  and  Kr  will  be  dropped  from  the 
notation  for  the  remainder  of  this  paper. 

Exactly  how  accurate  eqs  (1-3)  axe  depends  upon  the  particular  type  of  photodiode  and  the 
range  covered  by  the  parameters.  Figure  1 illustrates  the  accuracy  of  eq  (1)  for  the  1337 
type  photodiode  (modeled  as  described  in  Paxt  II  of  this  paper)  over  the  400  to  900  nm 
spectral  region  for  an  internal  quantum-deficiency  spectrum  that  is  a practical  upper  limit 
for  this  type  of  photodiode.  Figure  1 compares  the  internal  quantum-deficiency  spectra 
(^(A,  0,  71264  cm/s,  oo)  and  ^(A,  0, 0, 1 ms)  with  their  sum,  and  with  the  internal  quantum- 
deficiency  spectrum  ^(A,  0,71264  cm/s,  1 ms).  The  latter  spectrum  and  the  sum  of  the 
two  first  spectra  differ  by  less  than  2 x 10”®  over  the  entire  400  nm  spectral  region,  and 
cannot  be  distinguished  at  the  scale  of  the  figure.  (Any  errors  with  absolute  values  less 
than  2 x 10“®  are  negligible  for  the  purposes  of  this  paper.) 

Figure  2 illustrates  the  accuracy  of  eqs  (1-3)  over  the  400  to  900  nm  spectral  region  for  the 
case  where  r,.  = oo.  The  spectra  in  figure  2 were  calculated  for  the  range  of  parameters 
listed  in  table  1.  Cases  1-3  in  that  table  are  based  on  the  front-region  doping  distribution 
iV’^(a;)/iV’^(0)  shown  as  the  dashed  curve  in  Figure  3,  and  Case  4 is  based  on  that  shown 
as  the  solid  curve  in  the  same  figure.  Of  six  doping  distributions  measured  for  different 
1337  type  photodiodes,  the  two  shown  in  Figure  3 produce  the  largest  change  in  shape  of 
the  internal  quantum-deficiency  spectrum. 

The  shapes  of  the  front-region  internal  quantum- deficiency  spectra  in  figure  2 are  compared 
in  figure  4 by  normalizing  all  of  the  spectra  to  0.01  at  440  nm  and  by  plotting  their 
differences  relative  to  the  normalized  spectrum  for  Case  2.  Cases  1 and  3 compared  with 
Case  2 shows  how  the  shape  of  Sf(X)  chtinges  when  its  magnitude  is  increased  by  over  a 
factor  of  three  in  association  with  either  a substantial  decrease  or  a substantial  increase  in 
Ng,.  (Cases  1 and  3 bracket  a number  of  cases  that  were  examined  in  which  Ng,  and  5 
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Figure  1.  The  quantum- deficiency  spectra  ^(A,  0,  71264  cm/s,  oo)  (dashed  line)  and  (5(A,  0.  0. 1 
ms)  (dot-dashed  Hne)  and  their  sum  (solid  line). 
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Figure  2.  The  quantum- deficiency  spectra  ^(A,  iV,,,  5,  r^)  for  the  conditions  (Case  1 - 
upper  dashed  line;  Case  2 - lower  dashed  line;  Case  3 - upper  solid  line;  Case  4 - lower 
solid  line)  listed  in  table  1. 
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Table  1.  Recombination-related  parameters  used  with  the  doping  distributions  shown  in 
figure  3 to  model  the  quantum-deficiency  spectra  of  Hamamatsu  1337  type  photodiodes. 


CASE 

N.. 

(cm-^) 

s 

(10^  cm/s) 

Tr 

(s) 

Na{i) 

1 

0 

0.71264 

OO 

dashed  curve 

2 

-3  X 10'^ 

1.83500 

oo 

dashed  curve 

3 

-1  X 10*’ 

100.000 

OO 

dashed  curve 

4 

-4  X 10*’ 

1.56000 

oo 

solid  curve 
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CONCENTRATION  (cm 


Figure  3.  The  front-region  doping  distributions  used  to  calculate  the  front-region 

interned  quantum-deficiency  curves  in  figure  2. 
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Figure  4.  Differences  between  the  internal  quantum-deficiency  spectra  of  figure  2 for  Cases 
1 (dashed  line),  3 (dotted  line),  and  4 (solid  line)  when  normalized  to  0.01  at  440  nm,  and 
that  for  Case  2 when  normalized  to  0.01  at  440  nm. 
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were  independently  varied  within  the  range  of  values  shown  in  table  1.) 

Case  4 compared  with  Case  2 shows  how  the  shape  of  <5/(A)  changes  when  the  dopant 
distribution  is  changed  from  the  dashed  curve  to  the  solid  curve  in  figure  4.  This  effect  is 
much  larger  than  that  associated  with  varying  N,,  and  5. 

The  three  spectra  in  figure  4 characterize  the  uncertainties  to  be  associated  with  the  diode- 
to-diode  variations  in  IV,,,  S,  and  Na{x)  when  any  one  of  the  spectra  in  that  figure  is  used 
as  <^/(A)  with  Kf  = 0.01  in  eqs  (1)  and  (3).  These  spectra  should  be  multiplied  hy  K f /O.Ol 
if  Kf  0.01.  For  the  purposes  of  this  paper,  the  (?(A,  iV,,,  5,  oo)  spectrum  corresponding 
to  Case  4 in  table  1 will  be  used  for  Sf[\)  in  eqs  (1)  and  (3). 

The  front-region  internal  quantum  deficiency  also  depends  upon  two  front-region  recom- 
bination mechanisms  not  considered  above.  These  are  the  Auger  and  band-to-band  mech- 
anisms that  become  important  in  the  heavily  doped  front  region  of  the  photodiode.  Both 
effects  were  included  in  all  of  the  spectra  described  so  far  by  using  the  default  cross  sections 
that  are  built  into  PC-lD  for  these  recombination  mechanisms.  The  PC-lD  default  values 
for  the  Auger  recombination  cross  sections  give  an  internal  quantum- deficiency  spectrum 
with  a value  of  9 x 10”®  at  400  nm.  The  shape  of  this  spectrum  is  well  modeled  by  ^/(A) 
[5].  In  fact,  an  uncertainty  of  ±100%  can  be  tolerated  in  the  Auger  cross  sections  without 
causing  sm  error  as  large  as  2 x 10“®  anywhere  within  the  400  to  900  nm  spectral  region. 

For  the  doping  profile  shown  as  the  solid  curve  in  figure  3,  the  PC-lD  default  values  for 
the  band-to-band  recombination  cross  sections  give  the  lifetimes  shown  in  figure  5 as  a 
function  of  distance  from  the  oxide-silicon  interface  in  the  photodiode.  These  lifetimes 
give  an  internal  quantum- deficiency  spectrum  with  a value  of  6 x 10“®  at  400  nm.  This 
spectrum  is  not  as  well  modeled  by  Sf{X)  as  is  the  spectrum  for  Auger  recombination,  but 
an  uncertainty  of  ±50%  in  the  band-to-band  recombination  cross  sections  can  be  tolerated 
without  raising  the  error  associated  with  this  uncertainty  to  2 x 10“®  anywhere  in  the  400 
to  900  nm  spectral  region. 

Now  consider  <^r(A).  Figures  6 and  7 illustrate  the  accuracy  of  eqs  (1)  and  (3)  over  the 
400  to  900  nm  spectral  region  for  the  case  where  5 = 0.  (Note  that  when  5 = 0,  the 
internal  quantum- deficiency  spectrum  is  independent  of  iV,,.)  The  spectra  in  figure  6 were 
calculated  for  = 1 ms  and  10  ms,  with  iV,,  z=  5 = 0.  For  the  purposes  of  this  paper, 
<^r(A)  will  be  set  to  zero  in  eqs  (2)  and  (3)  for  A < 610  nm,  and  the  spectrum  for  = 1 
ms  will  be  used  for  A > 610  nm. 

The  shapes  of  the  spectra  in  figure  6 are  compared  in  figure  7 by  normalizing  them  to 
0.002  at  860  nm  and  plotting  their  differences  relative  to  the  spectrum  for  = 1 ms. 
The  normalization  value  of  0.002  was  chosen  as  a practical  upper  limit.  The  difference 
spectrum  in  figure  7 characterizes  the  uncertainties  to  be  associated  with  the  variations  in 
rear-region  lifetime  when  the  spectrum  for  = 1 ms  is  used  as  with  Kr  = 0.002  in 

eqs  (2)  and  (3).  This  spectrum  should  be  multiplied  by  Kr /0.002  if  Kr  ^ 0.002. 

A convenient  function  of  wavelength  has  been  fitted  to  the  spectrum  corresponding  to  Case 
4 in  table  1,  and  another  to  the  spectrum  for  = 1 ms  in  figure  6.  These  functions  are 
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Figure  5.  The  default,  band-to-band  recombination  lifetime  profile  calculated  by  PC- ID 
for  the  doping  distribution  shown  as  the  solid  curve  in  figure  3. 
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INTERNAL  QUANTUM  DEFICIENCY 


400  500  600  700  800  900 


WAVELENGTH  (nm) 

Figure  6.  The  internal  quantum-deficiency  spectra  ^(A,  iV,,,  5,  r^)  for  Tr  = 1 ms  (solid 
line),  and  = 10  ms  (dashed  line)  with  Ng,  = 5 = 0. 
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Figure  7.  Differences  between  the  internal  quantum-deficiency  spectra  of  figure  6 when 
normalized  to  0.002  at  860  nm. 
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^,(A)  = 4o[l/.Y(A)-A2/r(A)], 


(5) 


respectively,  where 


Ar(A)  = exp(-A/Ao), 


(6) 


and 


y(A)=exp(-(A/A2f).  (7) 

The  appropriate  values  for  and  A2  for  eqs  (4)  and  (5)  are  listed  in  table  2. 

The  function  in  eq  (5)  becomes  negative  for  A < 620  nm,  in  which  case  eq  (5)  is  replaced 
with  ^r(A)  = 0 as  mentioned  earlier.  The  values  of  Aq  in  table  2 normalize  <^/(A)  to  0.01 
at  440  nm,  and  ^r(A)  to  0.002  at  860  nm,  respectively. 

With  Kf  = 0.5516405,  eq  (4)  fits  the  spectrum  for  Ng,  = —3  x 10^^  cm/s  in  figure  4 with 
a residual  standard  deviation  of  2.9  x 10“®  over  the  420  to  900  nm  spectral  region  and 
never  differs  from  that  spectrum  by  more  than  7 x 10~®  over  the  same  spectral  region. 
This  translates  to  less  than  10  ppm  when  the  internal  quantum  deficiency  is  normalized 
to  0.01  at  440  nm;  thus  the  error  associated  with  the  use  of  eq  (4)  is  negligible  for  the 
purposes  of  this  paper.  The  larger  variations  below  420  nm  are  negligible  with  respect  to 
uncertainties  already  identified  and  others  discussed  below. 

With  Kr  = 1.685081,  eq  (5)  fits  the  spectrum  shown  as  a solid  line  in  figure  6 with  a 
residual  standard  deviation  of  1.4  x 10~®  over  the  610  to  900  nm  spectral  region,  and  never 
differs  from  that  spectrum  by  more  than  2.5  x 10“®  over  that  same  spectral  range.  This 
translates  to  less  than  1.7  x 10“®  when  the  quantum  deficiency  is  normalized  to  0.002  at 
440  nm;  thus  the  error  associated  with  the  use  of  eq  (5)  is  also  negligible  for  the  purposes 
of  this  paper. 

3.  Errors  Associated  with  Quantum- Deficiency  Calculations 

Two  types  of  error  associated  with  the  quantum-deficiency  values  calculated  by  PC- ID  can 
be  distinguished:  1)  numerical  errors  associated  with  the  algorithms  used  by  PC- ID  and 
the  limited  number  of  finite  elements  available  to  PC-lD,  and  2)  physical  errors  associated 
with  approximations  and  simplifications  in  the  physical  models  used  with,  or  built  into 
PC-lD.  Both  types  are  considered  in  this  section. 

It  was  shown  in  Ref.  [5]  that  Sf{\)  is  well  approximated  by 


Table  2,  Values  of  the  parameters  to  be  used  in  eq  (4)  for  extrapolating  a short  wavelength 
internal  quantum-deficiency  value  to  longer  wavelengths  and  in  eq  (5)  for  extrapolating  a 
long  wavelength  internal  quantum-deficiency  value  to  shorter  wavelengths. 


Parameter 

eq  (4) 

eq  (5) 

^0 

1.38002  X 10^ 

2.02157  X 10-® 

Ai 

1.47530  X 10^ 

3.51764  X 105- 

A2 

8.07476  X 10-5 

0.0 

^0 

4.39312  X 10^  nm 

7.53991  X 10^  nm 

A2 

4.27998  X 10^  nm 

1.66889  X 10®  nm 
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(8) 


Sf(\) 


r[J(0)-J(a(A))] 

l-r7(a(A)) 


/:(2;)  exp(  — a(A)x)ff2;, 


(9) 


k{x)  = k{x)/ k{0)^ 


(10) 


k{x)  = l/{mo{x)Dm{x)), 


(11) 


T = S/Dm{0),  (12) 

where  q(A)  is  the  absorption-coefficient  spectrum  of  silicon,  x is  the  distance  from  the 
oxide-silicon  interface  toward  the  rear  of  the  photodiode,  7710(1)  is  the  equilibrium  minority- 
carrier  concentration  at  the  point  x in  the  front  region  of  the  photodiode,  and  Dm{x)  is 
the  equilibrium  diffusion  constant  for  the  minority  carriers  at  the  point  x.  Note  that  the 
sign  in  the  numerator  in  eq  (8)  is  correct  and  that  the  sign  is  incorrect  in  Ref.  ^5;.  Also, 
note  that  /(A)  becomes  very  insensitive  to  xi  as  xi  increases  beyond  some  critical  value. 
The  point  where  the  electron  and  hole  concentrations  cross  is  a convenient  choice  for  xi, 
but  smaller  values  would  work  as  well. 

From  the  point  of  view  of  PC-lD,  the  kernel  k{x)  in  eq  (11)  is  given  by 


fc(i)  = qMo{x)l{n'i,{x)kTnm{x)),  (13) 

where  Mo{x)  is  the  equilibrium  majority  carrier  concentration  in  the  front  region,  nt^{x) 
is  the  effective  intrinsic  carrier  concentration  [6],  which  in  PC-lD  includes  the  effect  of 
Fermi-Dirac  statistics,  is  the  minority-carrier  mobility  as  a function  of  position  in 

the  front  region,  k is  the  Boltzmann  constant,  and  T is  the  temperature  of  the  photodiode. 

Figure  8 compares  exp(— a(400  nm)  x)  with  the  equilibrium  functions  iVfo(a:),  7iic(i),  and 
as  calculated  by  PC-lD  over  the  range  0 < x < 200  nm  for  the  same  photodiode 
model  used  to  calculate  the  quantum-deficiency  spectra  in  figure  3 for  Case  2 of  table  1. 
Equations  (8-12),  figure  8,  and  figure  3 show  why  the  shape  of  Sf{X)  is  very  insensitive  to 
Ntg  while  being  quite  sensitive  to  the  shape  of  Na{x)/^a{0).  Since  the  width  of  the  accu- 
mulation layer  created  by  N,,  is  small  compared  to  the  distance  over  which  exp(  — a(A)x) 
changes  significantly  for  A > 400  nm,  changes  in  Ng,  have  the  effect  of  multiplying  the 
integral  in  eq  (9)  by  a constant  factor  independent  of  A > 400  nm.  On  the  other  hand,  the 
changes  in  Na(x)/Na{0)  shown  in  figure  3 extend  well  into  the  region  where  exp(  — q(A)x) 
varies  significantly  from  unity  for  A = 400  nm.  As  A is  increased  above  400  nm,  the  shape 
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Figure  8.  Comparison  of  Mq{x)  (dotted  line),  nie{^)  (dashed  line),  /im(aj)  (dot-dashed 
line),  and  exp(  — a(400  nm)  x)  (solid  line)  for  a 1337  type  photodiode  with  iV,,,5,  and 
set  as  defined  in  Case  2 of  table  1.  The  tildes  in  the  preceding  expressions  indicate  that 
these  expressions  have  been  normalized  to  unity  at  z = 0,  as  shown  for  k{x)  in  eq  (10). 
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of  Sf(X)  will  become  less  dependent  on  the  shape  of  iV^(a:)/iV^(0)  because  q(A)  decreases 
with  increasing  A, 

The  numerical  accuracy  with  which  PC-lD  calculates  Mo{x)  within  the  framework  of  the 
drift-diffusion  approximations  for  the  conditions  of  charge  accumulation  at  an  oxide-sihcon 
interface  as  shown  in  figure  8 has  been  studied  previously  and  found  to  be  quite  good  '7]. 
With  a uniform  doping  density  of  10^®  cm“^  and  an  oxide  charge  density  of  10^^  cm“^, 
PC-lD  overestimates  Mo{x)  by  about  9%  at  x = 0,  decreasing  to  within  ±1%  for  z > 1 
nm.  The  uncertainty  in  k{x)  due  to  this  source  of  error  will  be  modeled  by  replacing  k(x) 

by 


^1(2)  = fc(a;)/(l  - O.lz/xo), 


(14) 


for  X < xq,  and 


ki{x)  = A;(a;)/0.9,  (15) 

for  I > lo  in  the  kernel  of  the  integrzd  in  eq  (9),  where  xq  = 1 nm. 

Another  source  of  error  that  must  be  considered  is  the  numerical  accuracy  with  which 
PC-lD  wordd  calculate  the  steady-state  photocurrent  if  the  equilibrium  carrier  concen- 
tration were  exact.  The  uncertainty  associated  with  this  source  of  error  was  estimated 
as  follows:  1)  The  internal  quantum-deficiency  spectrum  was  cedculated  directly  from  the 
PC-lD  solution  of  the  steady-state  drift- diffusion  equations  for  the  total  current  flowing 
in  the  photodiode  for  the  conditions  described  by  Case  2 of  table  1.  2)  Equations  (8)- 
(13)  were  used  to  calculate  the  internal  quantum-deficiency  spectrum  from  the  equilibrium 
values  of  MQ{x)^nic(x),  and  /Xm(®)  calculated  by  PC-lD  for  the  same  conditions  as  in  1) 
above.  The  integration  in  eq  (9)  was  carried  out  using  a generalization  of  Simpson’s  rule 
for  nonequi distant  points.  3)  Both  spectra  were  normalized  to  0.01  at  440  nm,  and  the 
differences  calculated.  There  is  no  reason  to  believe  that  these  differences,  which  exceed 
2 X 10“^  only  at  400  nm  and  at  900  nm,  are  caused  by  errors  associated  with  PC-lD.  It  is 
just  as  likely  that  they  are  associated  with  eqs  (8-12)  since  these  equations  were  derived 
from  a number  of  assumptions  that  are  not  rigorously  satisfied.  The  important  point  is 
that  it  is  unlikely  that  either  calcidation  is  in  error  by  much  more  than  the  differences 
between  them.  This  verifies  the  numerical  accuracy  of  the  photocurrents  calculated  by 
PC-lD  to  the  level  of  accuracy  required  in  this  paper. 

There  are  a number  of  other  errors,  which  are  associated  with  incomplete  or  approximate 
physical  models  built  into  PC-lD,  that  also  affect  the  accuracy  with  which  PC-lD  cal- 
culates the  functions  defining  k{x)  in  eq  (13).  These  include  quantum-mechanical  effects 
(tunneling)  [8]  that  chajige  the  shape  of  Mo{x)  from  that  calciilated  from  the  drift -diffusion 
approximation  for  x < 1 nm,  band-gap  narrowing  effects  [9]  associated  with  heavy  doping 
and  large  surface-fields  [10]  that  cause  nie{x)  [6]  to  depend  upon  position  in  the  photodi- 
ode, and  carrier-carrier  and  carrier-ion  scattering  mechanisms  that  cause  to  depend 
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upon  position  11^.  It  is  beyond  the  scope  of  this  paper  to  investigate  each  of  these  sources 
of  error  in  detail,  but  the  uncertainty  in  k[x)  due  to  this  source  of  error  will  be  modeled 
by  replacing  k{x)  by 


k2[x)  = k{x){l  - O.llog{k{x)))  (16) 

in  the  kernel  of  the  integral  in  eq  (9).  This  allows  an  error  of  10%  per  decade  change  in 
k{x). 

The  differences  between  the  internal  quantum-efficiency  spectrum  calculated  from  eqs  (3- 
13)  by  replacing  k{x)  with  either  ki{x)  or  k2{x)  in  eq  (9)  and  the  spectrum  calculated  using 
k(x)  were  calculated  after  all  three  spectra  were  normedized  to  0.01  at  440  nm.  Neither 
difference  exceeded  2 x 10“®  over  the  entire  400  to  900  nm  spectral  region.  Therefore,  ail 
of  these  sources  of  error  are  considered  negligible  for  the  purposes  of  this  paper. 

Another  source  of  error  that  is  obvious  in  eqs  (8-12)  is  the  uncertainty  associated  with  the 
absorption-coefficient  data  used  in  the  simulations.  As  described  in  Part  I,  the  absorption- 
coefficient  data  were  calculated  from  an  equation  [12]  fitted  to  the  data  of  Weakliem  and 
Redfield  [13].  An  alternate  set  of  absorption-coefficient  data  was  described  by  Philipp  14  . 
Figure  9 plots  the  differences  between  the  internal  quantum-deficiency  spectra  normalized 
to  0.01  at  440  nm  for  Case  2 of  table  1 when  based  on  the  absorption-coefficient  data  in 
Ref.  [14]  and  when  based  on  the  equation  in  Ref.  [12].  These  differences  will  be  used  as  the 
estimated  uncertainties  arising  from  the  uncertainties  in  the  silicon  absorption-coefficient 
spectrum  when  Sf{X)  is  used  to  extrapolate  an  internal  quantum  deficiency  of  0.01  from 
440  nm  to  longer  wavelengths  with  a 1337  type  photodiode.  A similar  calculation  was 
carried  out  for  the  case  where  5 = 0 and  = 1 ms.  Both  spectra  were  normadized  to 
0.002  at  860  nm.  The  differences,  which  are  shown  in  figure  10,  are  taken  as  the  estimated 
uncertainties  arising  from  the  uncertainties  in  the  silicon  absorption-coefficient  spectrum 
where  Sr{X)  is  used  to  extrapolate  an  internal  quantum  deficiency  of  0.002  from  860  nm 
to  shorter  wavelengths  with  a 1337  type  photodiode. 

Figures  11  and  12  plot  the  quadrature  sum  of  the  difference  spectra  plotted  in  figures  4 and 
9 and  in  figures  7 and  10,  respectively.  The  spectrum  in  figure  11  plots  the  estimated  (one 
standard  deviation)  uncertainty  as  a function  of  wavelength  that  is  associated  wdth  the 
use  of  <^/(A)  in  eq  (4)  to  extrapolate  a measured  internal  quantum  deficiency  of  0.01  from 
440  nm  to  longer  wavelengths.  The  major  sources  of  error  contributing  to  the  uncertainty 
spectrum  in  figure  11  are  the  diode-to-diode  variations  in  NaM/^aW  and  the  uncertainty 
in  the  absorption-coefficient  spectrum  for  silicon,  but  a number  of  other  sources  of  error 
were  identified  with  the  help  of  eqs  (8-12),  and  any  that  produced  differences  greater  than 
2 X 10“®  for  400  nm  < A < 900  nm  axe  included  in  figure  11. 

The  spectrum  in  figure  12  plots  the  estimated  (one  standard  deviation)  uncertainty  as  a 
function  of  wavelength  that  is  associated  with  the  use  of  ^r(A)  in  eq  (5)  to  extrapolate  a 
measured  internal  quantum  deficiency  of  0.002  from  860  nm  to  shorter  wavelengths.  The 
major  sources  of  error  contributing  to  the  uncertainty  spectrum  in  figure  12  are  the  change 
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Figure  9.  Differences  between  the  internal  quantum- deficiency  spectrum  for  a 1337  type 
photodiode  satisfying  the  conditions  defined  in  table  1 for  the  absorption-coefficient  data 
of  Ref.  [14]  and  for  the  absorption-coefficient  data  of  Ref.  [12]. 
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Figure  10.  Differences  between  the  internal  quantum-deficiency  spectra  for  a 1337  type 
photodiode  with  5 = 0 and  = 1 ms  for  the  absorption-coefficient  data  of  Ref.  [14:  and 
for  the  absorption-coefficient  data  of  Ref.  [12]. 
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Figure  11.  Sum  in  quadrature  of  the  curves  in  figures  4 and  9.  This  curve  is  the  nominal 
uncertainty  associated  with  the  use  of  eq  (4)  in  eq  (2)  to  extrapolate  a quantum  deficiency 
of  0.01  at  440  nm  to  any  other  wavelength  between  400  and  900  nm.  Notice  that  the  sum 
is  plotted  on  different  scales  below  and  above  440  nm. 
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Figure  12.  Sum  in  quadrature  of  the  curves  in  figures  7 and  10.  This  curve  is  the  nominal 
uncertainty  associated  with  the  use  of  eq  (5)  in  eq  (2)  to  extrapolate  a quantum  deficiency 
of  0.02  at  860  nm  to  any  other  wavelength  between  400  and  900  nm.  Notice  that  the  sum 
is  plotted  on  different  scales  below  and  above  440  nm. 
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in  shape  of  (5r.(A)  with  and  the  uncertainty  in  the  absorption-coefhcient  spectrum  for 
silicon.  All  other  sources  of  error  considered  were  negligible  with  respect  to  2 x 10“^  for 
400  nm  < A < 900  nm. 

No  equivalents  to  eqs  (8-12)  were  used  to  guide  the  error  analysis  for  <5r(A).  The  results 
of  the  uncertainty  analysis  of  <^/(A)  show  that  numerical  accuracy  of  PC-lD  will  be  a 
negligible  source  of  error.  This  leaves  the  errors  associated  with  the  approximations  in  the 
physical  models  that  are  built  into  PC-lD.  Some  of  these  were  tested,  but  no  tests  were 
possible  for  the  use  of  a single  SRH  trap  level  to  model  the  recombination  in  the  rear  region 
of  the  photodiode.  This  approximation  could  be  a non-negligible  source  of  error;  it  was 
shown  in  Part  I that  this  approximation  introduced  errors  of  the  order  of  the  effect  being 
modeled  in  nonlinearity  simulations  on  a UV444B  type  photodiode.  A more  complete 
model  might  make  depend  upon  position  in  the  photodiode,  which  might  modify  the 
shape  of  <^r(A)  and  add  an  uncertainty  that  is  not  included  in  figure  12. 

4.  Comparison  with  Experiment 

To  use  the  results  developed  thus  far,  it  is  necessary  to  define  the  internal  quantum- 
deficiency  interpolating  function  ^(A),  and  to  assume  that  it  exactly  satisfies  eqs  (1-3), 
which  can  be  rewritten  in  more  compact  form  as 


«(A)  = S.(Af) 


^/(A) 


+ ^i(Ar) 


#r(A) 

5,(A.)J  i.(A.)’ 


(17) 


where  <^*(A/)  and  ^x(Ar)  are  the  measured  values  of  the  internal  quantum  deficiency  at  the 
wavelengths  Af  and  A^.,  respectively;  S/(A)  and  ^r(A)  are  defined  in  eqs  (4)  and  (5),  respec- 
tively, and  the  constajits  used  in  these  equations  are  defined  in  table  2.  The  uncertainty 
associated  with  ^(A)  due  to  the  sources  of  error  considered  in  the  last  two  sections  of  this 
paper  can  be  obtained  by  summing  in  quadrature  each  of  the  terms  on  the  right-hand  side 
of  the  differential  of  eq  (17).  That  differential  is  given  by 


dS(A)  = dS^(Af) 


' SfW ' 


-|-  d6x{Ar) 


+ dSf(A) 


-b  dSr{A) 


«r(A)  1 

L^r(Ar)J 

L «r(Ar)J 

^*(A/)] 

«/(A/)J 

«.(Ar)l 

f, 

«r(A.)J 

L «r(Ar)J 

(18) 
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The  differentials  dSf{X)  and  dSr{X)  are  plotted  in  figures  11  and  12  for  Xf  = 440  nm,  and 
for  = 860  nm,  respectively. 

There  are  not  many  data  available  against  which  to  test  eqs  (17)  and  (18),  but  there  are 
some.  For  instance,  Zalewski  and  Hoyt  [4]  have  reported  the  internal  quantum  efficiencies 
at  441.6  nm  and  the  spectral  responsivities  and  absorptances  at  633.0  nm  for  five  multiple 
reflection  (trap)  [4,15-17]  radiometers  based  on  Hamamatsu  1337  photodiodes  at  441.6 
and  633.0  nm.  Similarly,  Fox  [17]  has  reported  the  average  internal  quantum  efficiency 
of  ten  trap  radiometers  also  based  on  Hamamatsu  1337  photodiodes  at  six  Ar^  ion  laser 
lines,  and  he  has  also  measured  the  oxide-bias  correction  for  a single  Hamamatsu  1337 
photodiode  at  four  Ar'^  ion  laser  lines. 

Zalewski  and  Hoyt  [4]  presented  their  data,  which  are  summarized  in  bold  face  type  in  table 
3,  as  an  intercomparison  of  spectral  responsivity.  Since  the  internal  quantum  deficiency 
is  the  quantity  of  interest  in  this  paper,  it  is  necessary  to  compute  this  quantity  from  the 
data  reported  by  Zalewski  and  Hoyt  as  shown  in  table  3.  The  corrections  for  nonlinearity 
that  Zalewski  and  Hoyt  derived  from  reverse-bias  measurements  also  eliminate  any  effects 
of  recombination  in  the  rear  region  of  the  photodiodes. 

Fox’s  [17]  oxide-bias  data,  which  are  summarized  in  table  4,  were  measured  at  10  V,  and 
were  reported  as  a fractional  increase  in  photocurrent  I given  by  7o  = J(10F)//(0)  — 1. 
The  oxide-bias  experiment  reported  in  Part  II  shows  that  the  internal  quantum  efficiency 
at  zero  bias  8^  is  given  by 


6,  = 1 -1/(1  +1.137o), 


(19) 


where  the  maximum  oxide  bias-voltage  is  10  V.  The  average  internal  quantum  efficiencies 
reported  by  Fox  are  also  shown  in  table  4,  and  the  internal  quantum  deficiencies  calculated 
from  them  are  also  shown  there. 

Equation  (17)  was  normalized  to  the  results  of  Zalewski  and  Hoyt  at  A/  = 441.6  nm,  and 
the  predictions  of  that  equation  are  compared  with  the  results  of  Zalewski  and  Hoyt  at  633 
nm  in  table  5.  Note  that  ^*(At.)  = 0 for  this  data  set  due  to  the  application  of  reverse  bias. 
Equation  (17)  was  also  normalized  to  the  oxide-bias  results  of  Fox  at  A/  = 468.18  nm,  and 
the  predictions  of  that  equation  are  compared  with  the  oxide-bias  results  in  table  6.  Again 
note  that  <^*(At.)  = 0 because  oxide-bias  measurements  are  not  sensitive  to  recombination 
in  the  rear  of  the  photodiode. 

The  uncertainty  estimates  for  eq  (17)  in  tables  5 and  6 made  use  of  the  data  in  figure  11 
even  though  those  data  apply  to  A/  = 440  nm  rather  than  441.6  or  468.2  nm.  This  does  not 
distort  the  estimated  uncertsdnty  significantly  because  the  terms  are  added  in  quadrature, 
and  the  first  term  on  the  right-hand  side  of  eq  (18)  is  larger  than  the  second  term  for 
both  sets  of  data.  Fox’s  data  point  at  406.85  nm  was  included  in  table  6 even  though  the 
uncertainties  associated  with  extrapolating  8f(X)  to  wavelengths  shorter  than  Xf  are  very 
unfavorable  for  high-accuracy  applications.  The  good  agreement  at  this  wavelength  must 
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Table  3.  Average  values  of,  and  uncertainty  estimates  for  the  internal  quantum  efhciency 
e/,  spectral  responsivity  R,  one  minus  the  reflectance  (1  — p)  and  nonlinearity  correction 
(l-NL)  reported  in  Ref.  [5]  (bold  face  type)  at  wavelength  A for  five  multiple-reflection 
radiometers  based  on  1337  type  photodiodes,  and  internal  quantum-deficiency  values 
calculated  from  the  reported  data. 


A 

(nm) 

R 

(A/W) 

(1~P) 

(l-NL) 

8,{\) 

441.6 

+0.99347 

±0.00040 

±0.00653 

±0.00040 

633.0 

+0.50840 

+0.99704 

±0.00009 

+0.99972 

±0.00014 

±0.00097 

±0.00004 

±0.00023 

Table  4.  Oxide-bias  data  7^  reported  in 
average  internal  quantum-efficiency  data 

Ref.  [17]  for  a single  1337  type  photodiode,  the 
e*  reported  in  Ref.  [17]  for  ten  multiple-reflection 

radiometers  based  on  1337  type  photodiodes,  and  the  quantum  deficiencies  <^x(A)  calculated 
from  these  sets  of  data. 

A — oxide  bias  — — iqe  — 


(nm) 

7o 

^x(A) 

^x(A) 

406.85 

±0.0126? 

±0.00015 

±0.01404 

±0.00017 

±0.9867 

±0.0002 

±0.0133 

±0.0002 

468.18 

±0.0044? 

±0.00015 

±0.00495 

±0.00017 

±0.9934 

±0.0003 

±0.0066 

±0.0003 

568.35 

±0.9977 

±0.00020 

±0.0023 

±0.0002 

647.30 

±0.0007? 

±0.00015 

±0.00079 

±0.00017 

±0.9986 

±0.0002 

±0.0014 

±0.0002 

799.54 

±0.9989 

±0.0002 

±0.0011 

±0.0002 

859.07 

±0.0000? 

±0.00015 

±0.00000 

±0.00017 

±0.9987 

±0.0002 

±0.0013 

±0.0002 
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Table  5.  Comparison  of  predictions  of  eq  (17) 

with  measured  internal 

quantum-deficiency 

data  <^z(A)  in  table  3 for  Xf 

= 441.6  nm  and  <^z(Ar)  = 0. 

Prediction 

Measured 

Predicted 

Difference 

Wavelength 

Value 

Value 

(nm) 

Ref.  [4] 

eq  (17) 

633.0 

+0.00097 

+0.00089 

-0.00008 

±0.00023 

±0.00007 

^0.00024 

Table  6.  Comparison  of  predictions  of  eq  (17)  with  internal  quantum-deficiency  data  (^^(A) 
derived  from  oxide-bias  data  in  table  4 for  Xf  = 468.18  nm  and  ^j.(Ar)  = 0. 

Prediction 

Measured 

Predicted 

Difference 

Wavelength 

(nm) 

Value 

Ref.  [15] 

Value 
eq  (17) 

406.85 

+0.01404 

+0.01321 

-0.00083 

±0.00017 

±0.01144 

±0.01144 

647.30 

+0.00079 

+0.00090 

+0.00010 

±0.00017 

±0.00030 

±0.00034 

859.07 

+0.00000 

+0.00015 

+0.00015 

±0.00017 

±0.00005 

±0.00018 
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be  considered  fortuitous  considering  the  uncertainty  associated  with  the  predicted  value. 
The  differences  at  the  other  wavelengths  never  exceeds  0.00015,  and  they  fall  within  the 
estimated  uncertainty  for  the  differences. 

Table  7 and  figure  13  compare  the  internal  quantum-deficiency  results  of  Fox  !17]  with  the 
predictions  of  eq  (17)  when  normalized  at  A/  = 468.18  nm  and  at  A^  = 859.07  nm.  The 
data  in  figures  11  and  12  that  apply  for  Ay  = 440  nm  and  for  A^  = 860  nm,  respectively,  are 
used  with  eq  (18)  to  estimate  the  uncertainties  assigned  to  the  predictions  of  eq  (17).  The 
difference  at  406.85  nm  is  much  larger  than  that  obtained  with  the  oxide-bias  data  but  still 
falls  within  the  estimated  uncertainty.  Because  the  oxide-bias  experiment  is  not  sensitive 
to  the  quantum  yield  for  electron-hole  pair  production,  it  is  tempting  to  imagine  that  the 
discrepancy  between  the  oxide  bias  and  interned  quantum  efficiency  reflects  the  fact  that 
the  quantum  yield  for  the  1337  type  photodiode  is  greater  than  unity  at  406  nm.  However, 
the  fact  that  this  discrepancy  is  not  statistically  significant  at  the  one-sigma  level  shows 
that  this  conclusion  cannot  be  drawn  from  the  data  presented  here.  The  differences  at  the 
other  wavelengths  never  exceed  0.0003,  and  only  fall  outside  the  estimated  uncertainties 
for  the  differences  at  799.54  nm. 

The  results  shown  in  tables  5 through  7 and  figure  13  verify  that  the  internal  quantum 
deficiency  spectrum  of  multiple-reflection  radiometers  based  on  Hamamatsu  1337  photo- 
diodes c£Ln  be  interpolated  over  the  440  to  860  nm  spectral  region  from  only  two  measured 
quantum  deficiencies,  one  at  each  end  of  the  region.  This  greatly  reduces  the  number 
of  measurements  needed  for  high-accuracy  spectral  calibrations  of  these  photodiodes  and 
radiometers.  Furthermore,  the  theoretical  uncertainty  analysis  suggests  that  the  uncer- 
tainties in  the  interpolated  values  will  be  less  than  or  equal  to  the  uncertainties  in  the 
measured  values,  at  least  as  long  as  the  latter  are  greater  than  ±0.0001.  The  results  in 
tables  5 to  7 and  figure  13  do  not  contradict  this  idea,  but  the  experimental  data  are  not 
accurate  enough  to  confirm  it.  On  the  other  hand,  the  experimental  data  do  confirm  that 
an  uncertainty  of  ±0.0003  is  obtainable. 

5.  Conclusion 

Part  III  has  shown  that  Version  2 of  the  semiconductor-device  modeling  program  PC- ID 
can  be  used  to  model  the  spectral  shape  of  the  internal  quantum  deficiency  of  1337  type 
photodiodes  with  very  small  uncertainty  over  the  440  nm  spectral  region.  The  largest 
uncertainties  are  caused  by  uncertainties  in  the  absorption- coefficient  spectrum  of  silicon, 
diode-to-diode  variations  in  the  shape  of  the  front  region  dopant  distribution,  smd  diode- 
to-diode  variations  in  the  position-independent  rear-region  minority- carrier  lifetime.  It 
was  not  possible  to  estimate  the  uncertainty  caused  by  approximating  the  rear-region 
minority-carrier  lifetime  as  being  independent  of  position. 

Simple  formulae  were  derived  that  allow  high-accuracy  internal  quantum-efficiency  cali- 
brations on  1337  type  photodiodes  to  be  interpolated  from  measured  values  at  the  ends  of 
the  440  to  860  nm  spectral  region  over  the  interior  of  that  region.  The  uncertainties  that 
can  be  obtained  with  these  formulae  are  comparable  to  those  that  can  be  obtained  with 
the  highest  accuracy  measurements  currently  available. 
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Table  7.  Comparison  of  predictions  of  eq  (17)  with  measured  internal  quantum-deficiency 
data  in  table  4 for  A / = 468.18  nm  and  A^.  = 859.07  nm. 


Prediction 

Measured 

Predicted 

Difference 

Wavelength 

Value 

V^alue 

(nm) 

Ref.  [15] 

eq  (17) 

406.85  nm 

+0.0133? 

+0.01744 

4-0.00414 

±0.00025 

±0.01000 

±0.01000 

568.35  nm 

+0.0023? 

+0.00224 

-0.00006 

±0.00030 

±0.00007 

±0.00031 

647.30  nm 

+0.0014? 

+0.00121 

-0.00019 

±0.00020 

±0.00004 

±0.00020 

799.54  nm 

+0.0011? 

+0.00082 

-0.00028 

±0.00021 

±0.00009 

±0.00023 
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Figure  13.  Comparison  of  measured  data  of  Ref.  [17]  (open  circles)  with  eq  (17)  (solid 
line)  normalized  to  pass  through  the  measured  point  at  468.18  and  859.1  nm.  At  these 
wavelengths  the  error  bars  are  the  uncertainties  in  the  measured  values;  at  the  other 
wavelengths  the  error  bars  are  the  quadrature  sum  of  the  uncertainties  in  the  measured 
values  and  those  associated  with  the  use  of  eq  (17). 
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